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Chapter 1 



Introduction to quantum 
entanglement 



The superposition principle — one of the several postulates of quantum mechanics — produces conse- 
quences that deviate from the predictions of classical mechanics. Consider, e.g., a pair of spin-1/2 
particles, each of which can inhabit an up (f)) or down (],)) spin state relative to some spin quan- 
tization axis or any linear combination of |t) and |J,). Then two spin-1/2 particles can, e.g., be in a 
state tti) or |J,|), in which each spin has a definite direction. However, they are also allowed to be 
in the singlet state 



in which neither spin possesses a definite direction. States such as |*~) are said to be entangled, 
a term coined by Erwin Schrodinger in 1935. The two spins, which we shall call A and B, can be 
spatially far from each other, even — in principle — at galactic separations. 

Suppose that a Stern-Gerlach measurement is performed on A and the result |t) is obtained. 
Then the spin state of B is in a definite state [[). However, if the result [[) is obtained for A then 
B, as required by quantum mechanics, is in the state [t)- Quantum mechanics appears to imply a 
nonlocal correlation between the two spins. More specifically, if one consider observables for the 
state 1^*"), such as and a^, one obtains the following predictions from quantum mechanics: 



Therefore, quantum mechanics can have the consequence that the expectation value of the product 




(1.1) 



(ai®aS) = -l, (a^) = (a|)=0. 



(1.2) 
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of the two observables is not necessarily equal to the product of the expectation values of the two 
observables, e.g., 



Actually, the results in Eq. (|1.2p can be easily explained by a classical theory, e.g., the ensemble 
of the pairs of spins in the measurement is a statistical equal mixture of || and ||. The question of 
whether quantum mechanics can predict something that cannot be explained by classical theories 
remained unanswered until in 1964, when John S. Bell came up with an inequality that all classical 
local theories (usually called local hidden variable [LHV] theories; local in the sense that the 
measurement outcome on one side should not influence the other side) obey, whereas it can be 
violated by quantum mechanics. Clauser, Horne, Shimony, and Holt (CHSH) [Ij later derived an 
inequality based on correlations, such as that in Eq. (jl.2p . that all classical local theories must 
satisfy: 



where a and a' (b and b') are unit vectors, representing pairs of different orientations of the Stern- 
Gerlach apparatus at A (B). This Bell-CHSH inequality can be violated by the singlet state |^~), 



violation has been demonstrated experimentally, e.g., using photon polarization states [21 [3]. 

So far, we have been discussing the singlet state. What about other entangled states? Gisin [4] 
showed that any non-product two-spin-1/2 state indeed violates the CHSH inequality. This non- 
product, unfactorizable property of wavefunctions is entanglement, the characteristic trait of quan- 
tum mechanics, and is responsible for the deviation from classical theories. 

Although entanglement was initially associated with a rather philosophical debate over the 
foundations of quantum mechanics, it has recently been discovered to be a useful resource. For 
example, for quantum cryptography Ekert [5] found that violation of the CHSH inequality can be 
used as a test of security in the process of the random key distribution via a spin singlet state. 
Bennett and Wiesner [6] found that transmission of a two- level quantum state (i.e., qubit), which 
is initially maximally entangled with another two-level system at the receiving end, can encode 
classical information of two bits ( "super dense-coding" ) . Bennett and co-workers [7] further found 




(1.3) 



{a ■ a <Si a ■ b) + {a ■ a a ■ 9) + {a ■ a' a ■ b) - {a ■ a' <Si a ■ V) < 2, 



(1.4) 



with the maximum value of the left-hand side being 2\/2 (see Appendix |A] for more details). The 
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that with the shared maximal entanglement (i.e., a singlet state), an unknown two- level state can be 
faithfully reconstructed via communication of only two classical bits ("quantum teleportation"). All 
these tasks are made possible by entanglement. Although it is not yet clear whether entanglement 
is necessary for the speed-up of quantum algorithms, such as Shor's factoring [8j and Grover's 
searching [9j, it has been established that entanglement does enable quantum computation. This 
is well illustrated by the so-called one-way quantum computer, due to Rossendorf and Brigel |10j . 
in which an initial highly entangled state (specifically, a "cluster" state), together with subsequent 
local measurements alone, allows efficient execution of quantum computation. 

Although violation of Bell inequalities is a necessary and sufficient signature of entanglement 
in the pure-state setting, the situation is more subtle in the setting of mixed states. What do we 
mean by a mixed state? When is a mixed state entangled? Let us begin with a familiar setting 
for density matrices, quantum statistical mechanics, where a Hamiltonian 7i leads to the density 
matrix p, given by 



where Z = Tr(e~^ ), [3 = l/ksT, and \n) is the energy eigenstate with energy eigenvalue En- A 
more general description of a system than by wavefunctions is thus provided by a density matrix, 
which can be regarded as a probabilistic mixture of pure states, hence, the name mixed state: 



with < Pi < 1 and ^^Pi = 1. Both a wavefunction and a density matrix p will hereafter 
be simply referred to as a state. The density a matrix description is necessary when the system 
interacts with an environment or when we only have access to part of a larger system. We remark 
that the decomposition in Eq. (jl.6p is by no means unique. For example, the "white" distribution 
^1 can be expressed as 




(1.5) 



n 




(1.6) 



T)(TI + 2U)(i 



(1.7) 



or 
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where |±) = (|t) ± |i))/\/2. Further more, there is no requirement that the pure states in the 
decomposition be orthogonal to one another. 

The idea of a mixed entangled state of two or more parties is naturally extended from the pure- 
state state that allows no decomposition into a mixture of factorizable pure states. If a 
mixed state does have such a decomposition, it is said to be separable or unentangled. Despite its 
seemingly innocuous definition, the question of whether or not a mixed state is entangled turns out 
to present deep mathematical challenges. There have been proposed several useful criteria (usually 
called separability criteria) that can (but not always) determine whether or not a given state is 
entangled. It has been found that certain mixed states, although entangled, do not violate any Bell 
inequality and, even more surprisingly, allow a classical, local description [TT]. Thus we see that 
entanglement in the mixed-state scenario is a much richer and subtler phenomenon than it is for 
pure states. 

As Bell inequalities do not, in general, completely reveal entanglement for mixed entangled 
states, other approaches to quantifying entanglement have emerged. These include (i) entanglement 
as a quantifiable resource: entanglement of distillation and entanglement cost; (ii) information- 
theoretic considerations: relative entropy of entanglement and related measures; and (iii) other, 
mathematical approaches, including the central theme of this dissertation. Approach (i) is perhaps 
the most natural way to quantify entanglement. However, as we shall see later, it has so far been 
limited to the settings of two parties (i.e., bi-partite systems), and there are major difficulties in 
extending it to multi-paritite systems. Approaches (ii) and (iii) are thus seen as indispensible 
for providing a better understanding of entanglement in various settings, especially multi-partite 
ones. Although the measure that we shall focus on for most of this dissertation belongs to the 
mathematical approach, later we shall show that it is nevertheless related to other entanglement 
and physical properties. 

As one of the central themes in the entanglement theory is to quantity the degree of entangle- 
ment, in the following, we shall introduce several important, standard measures of entanglement. 
We shall conclude this chapter with an overview of the dissertation. 
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1.1 Separability and entanglement 

Let us define precisely whether a state is entangled or not. A state is entangled if it is not separable. 
A bi-partite state, describing parties A and B, is separable (or unentangled) if and only if 
can be expressed as 

i 

where the {pf^^^Ys are local states of A(B), which can be either pure or mixed, and {pij's are 
probabilities with < < 1 and 'Y^^Pi = 1- Such a sum, in which the weights are non-negative, is 
called a convex sum. The generalization to multi-partite states involves including more parties: 



pt''-'' = Y.p^pt®pf®---®pf- (1-10) 



In the present chapter we shall focus on bi-partite states. If a bi-partite state cannot be written as 
a convex sum of direct products of density matrices then it is entangled. However, this definition 
does not offer a practical way of determining separability or entanglement. 

Peres |12j proposed a very simple but useful criterion for separability. As a density matrix is 
Hermitian and positive semi-definite, its transpose is still a valid density matrix. If we take the 
transpose of the matrices {/cfj's in Eq. (jl.9p . the resulting matrix, denoted by pj^, still contains 
non- negative eigenvalues. The operation is called partial transpose and can be defined for any 
bi-partite state: 

P=Y. PirM\4 ® e]^>(efc ®ef\^ p^^=Y, Pihkj\ef )(e^ ef |, (1.11) 

i,j,k,l 

where |ef^ (gi e^) = |ef^) (Xi is the product basis used to represent the density matrix, and the 
underscores are used to highlight the changes under the partial transpose. Thus we have that if 
the state is separable, its partially transposed matrix has non-negative eigenvalues (usually called 
PPT). Said equivalently, if the state is not PPT under the partial transpose, the state must be 



entangleqj. This is the Peres positive partial transpose (PPT) criterion for separability [12 



^In general, if a state has PPT, no entanglement can be distilled out from it 19 . But the state can be either 
unentangled or entangled. When the state has PPT and is also entangled, it is called a bound entangled state. 
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Let us examine the example of a singlet state | ) . When written in the form of density matrix 
in the basis {|tT) , |Ti) , lit) , it corresponds to the density matrix 




i 




-i 











(1.12) 



\0 oj 



The partial transpose takes it to 








2 



i 



V- 



i 

-i 



(1.13) 



which has one negative eigenvalue, —1/2. Thus, via the PPT criterion we see that |^'~) is entangled. 

In general, this PPT criterion is necessary but not sufficient for establishing separability. How- 
ever, it was shown by Horodecki and co-workers [13j that PPT is sufficient in the cases of (8) 
(two-qubit) and (7^0 (qubit-qutrit) systems: if two-qubit or qubit-qutrit states obey PPT, they 
are separable. On the other hand, if PPT is violated, the state is entangled. The extent to which 
a state violates PPT is manifested in the negative eigenvalues of the partially transposed density 
matrix, and can be used as a measure (not just an identifier) of entanglement; this measure is called 
the negativity |14[|15| . Following Zyczkowski and co-workers [14j, we define the negativity M to be 
twice the absolute value of the sum of the negative eigenvalues: 



Af{p) = 2max(0, -Aneg), 



(1.14) 



where Aneg is the sum of the negative eigenvalues of p^^ and the factor of two is a normalization 
chosen such that the singlet state has Af = 1. 
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1.2 Entanglement of distillation 

The notion of the entanglement of distillation was introduced by Bennett and co-workers \16\ [T7] to 
give an operational definition of the degree of entanglement. Suppose p represents the state of two 
particles possessed by two parties (usually referred to as Alice and Bob) separated by some distance. 
A way to envisage the degree of entanglement that p has is to ask how useful p is compared to a 
standard state, such as any of the four Bell states: 



Here {|0) , |1)} represents an orthonormal basis of a two-level system, for instance, the z-component 
of the spin of a spin-1/2 particle, or the polarization of a photon. More specifically, given n copies 
of the state p shared between Alice and Bob, how many pairs, say k, of Bell states can be obtained 
if each of Alice and Bob is allowed to (i) perform any local operations (including measurement) 
on the particles he or she possesses and (ii) share with the other party classical information, e.g., 
the outcome of some measurements. These operations are called local operations and classical 
communication (LOCC). The asymptotic limit 



is called the entanglement of distillation [16\ I17j . In it, k is the average number of Bell states 
taken over different possibilities (due to measurement) of an optimal procedure. E^, quantifies the 
entanglement as a resource, using Bell states as a standard ruler. 

Let us illustrate the idea with an example. Suppose /? is a pure state corresponding to the ket 



with 9 € [0,7r/2]. For 9 = or 7r/2, it is not entangled. For the intermediate range of 9, the state 
is entangled, and maximally so at = 7r/4. Suppose that two copies of non-maximally entangled 




(1.15) 




(1.16) 




(1.17) 
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lipe) are shared between Alice and Bob: 



l^e) 12 ®\'<p0)u = (cos^|0i02) + sin^|lil2))®(cose|0304)+sine|l3l4)) 

= cos2^|0i020304) + V2cos^sme^(|0i02l3l4) + II1I2O3O4)) + Sin2e|lil2l3l4), 

v2 

where Ahce has particles 1 and 3, whereas Bob has 2 and 4. Their joint goal is to extract a Bell 
state under LOCC. 

Both parties can perform any local operations allowed by quantum mechanics. A possible 
operation is to measure the number of his/her particles in state |1) (e.g., the z-component of 
total angular momentum). If Alice measures the number of I's of her particles to be or 2 then 
the resulting state (IO1O2O3O4) or II1I2I3I4)) is unentangled. She needs to tell Bob to abort the 
operation, as there is now no entanglement to extract. But with probability 2cos^ 0sin^ 9 she gets 
the state 



^(|0l02l3l4) + |lll20304)), 



(1.18) 



which is evidently entangled. But how do they establish from this a Bell state, say, | 

This time, Alice proceeds to perform a unitary transformation U on her particles and contacts 
Bob (which is when the classical communication takes place) and asks him to perform the same 
unitary transformation on his particles. Suppose that the unitary transformation they agree to 
perform is (also known as a CNOT operation) 



U 



10 
10 
1 
10 



(1.19) 



in the basis of {|00), |01), |10), |11)}. In particular, for it we have 



C/13IO1I3) = IO1I3), UislUOs) = II1I3), i724|02l4) = IO2I4), U24\120a) = {hU) ■ 



(1.20) 
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Entangled sources 



Bob 



-IT 



Figure 1.1: Entanglement of distillation. 



Then the joint state after the transformations becomes 



2 
4 



U13U24 (^-^(lOiOslsU) + II1I2O3O4)) 

^(|0l02l3l4> + |lll2l3l4)) 
1 



7=2 



|0l02) + |lll2))®|l3l4). 



(1.21) 



As particles 3 and 4 are not entangled with 1 and 2, what now needs to be done is that Alice throws 
away her particle 3 and Bob throws away his particle 4. Finally, they have distilled one maximally 
entangled pair -^(|0i02) + II1I2)) out of two non-maximally entangled pairs. The probability P 
of success is 2cos^^sin^0, i.e., on average they can distill k/n = = cos^^sin^^ Bell pairs per 
initial pair. 

The above discussion involves Alice and Bob dealing with two pairs at a time. In fact, it can 
be extended to the case where they can manipulate n copies at a time |16j . The average number 
of Bell pairs per initial pair can be dervied to be 



^ n 1 " 

E{n) = -J2 = - E ^(^) log2(Cr), 



k=0 

2 a\n~ki 



k=0 



p{k) = {cos' eT~\s\n^ef CI, 



(1.22a) 
(1.22b) 



where = n\/[k\{n — A;)!]. As the number n of copies approaches infinity. 



lim E{n) ^Eb = S{pa) 



(1.23) 



where pA = TrB|V'e)(V'e|i and S{p) = — Trplog2/0 is the von Neumann entropy of p. In the case of 
{i/je), its entanglement of distillation is Ed = /i(cos^ 6), where h{x) = — x log2(a;) — (1 — x) log2(l — x), 



9 



0.2 0.4 0.6 0.8 1 

r 



Figure 1.2: Mixed state distillation for Werner state p\\f^{r). States above the dashed line (i.e., 
r > 1/3) can be distilled by the procedure. 

i.e., is the Shannon entropy. The result 

£;d = -Tr(/>Alog2(/9A)), (1.24) 

is valid for any bi-partite pure state. 

In Fig. 11.11 we show a distillation scheme slightly modified from the two-pair example. In this 
modified scheme, Alice and Bob both perform the CNOT operation before the measurement. This 
transforms the initial state \tpe)i2 ® {^^9)34 ^ follows: 

l^e)i2®lV'e>34 ^ (cos2^|OiO2)+sin2^|lil2))|O3O4) + ^cosesin04s(|OiO2) + |lil2»|l3l4>. (1.25) 

v2 

If Alice and/or Bob then measures the third and/or fourth qubit, respectively, and the outcome is 
|1), they immediately obtain a Bell state shared between particles 1 and 2. If the outcome is |0), 
they get a slightly less entangled state, which they can store for a second trial of distillation. What 
we mean by this is that two pairs of the states 

^ 'cos2 6'|00) +sin2 6'|ll)), (1.26) 



\/ cos'^ 9 + sin 
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although less entangled than the original pairs of Eq. (jl.l7p . is distillable. Thus, this modified 
scheme performs slightly better than the original two-pair scheme. 

For mixed entangled states, there are very few cases for which Ej) is known. No general optimal 
distillation procedure is known for generic states. But a similar set-up to the one shown in Fig. 11.11 
(except that the measurement is performed in a different basis) does provide a way (although not 
optimal) to distill very general two-qubit states. For example, Bennett and co-workers [T7] have 
shown that after one step of the mixed-state distillation procedure, two initial pairs of the state 
(which is usually called the Werner state) 

pw+(r) = r|^+)(^+| + i^l, (1.27) 



will be transformed into one pair with a new parameter (see Fig. 11.2 

, _ 2r(l + 2r) 
^ ~ 3(1 + r2) ■ 



(1.28) 



Note that the larger the parameter r is, the higher entanglement the Werner state possesses. If 
r' > r, i.e., when r > 1/3 or equivalently the fidelity F = {^^\ pw |^^) > 1/2, the entanglement is 
said to be increased. Horodecki and co-workers [18] further showed that any entangled two-qubit 
state can be transformed into a state p^^^{r) with r > 1/3, and hence can be distilled via the 
scheme of Bennett and co-workers (also known as the BBPSSW scheme) |17j . 

However, the procedure is not optimal for an arbitrary state p, and it is generally rather difficult 
to compute E-q{p). Nevertheless, if E-£){p) > we say that the state p is distillable. We remark 
that there is a connection between the PPT criterion and the distillability of a bi-partite state. 
Horodecki and co-workers [19] found that if a state has PPT then it cannot be distilled. But the 
converse is not generally true. In fact, there is still no simple criterion to determine whether or not 
a state is distillable. 
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1.3 Entanglement cost and entanglement of formation 

The distillation is a process for concentrating entanglement from a large number of pairs with less 
entanglement into a small number of pairs with more (and ultimately maximal) entanglement. On 
the other hand, we can consider the converse process, which is usually called dilution. Given k pairs 
of Bell states shared between Alice and Bob, how many pairs n of a given state p can be obtained 
by local operations (including adding unentangled particles) and classical communication? The 
goal is to maximize the number n of copies of the output state p. The optimal ratio defines the 
entanglement cost [20j : 

Ecip) = lim (k/n). (1.29) 

k—tca 

As with i?D, Eq is very difficult to calculate for general mixed states, and is only known for a very 
few special cases. 

However, for pure states such as the state \ipg) discussed previously, Eq = — Tr(/>^ log2(/9yi)) , 
which equals Ej). The optimal way to realize this dilution process for the pure state is to utilize 
two techniques: (i) quantum teleportation, which we have introduced at the beginning and which 
simply says that a Bell state shared between two parties can be used to tranfer an unknown qubit 
state with certainty, and (ii) quantum data compression [21j, which basically states that a large 
message consistive of say n qubits, with each qubit on average being described by a density matrix 
PA, can be compressed into a possibly smaller number k = nS{pA) < n of qubits; and one can 
faithfully recover the whole message, as long as n is large enough. For more detail of quantum data 
compression, see Appendix iBl 

With these two tools in hand, Alice can first prepare n copies of {ipg) (2n qubits in total) 
locally, compress the n qubits to k qubits that she will "send" to Bob, and teleport the compressed 
k qubits to Bob using the shared k Bell states. Bob then decompresses the k qubits back to the 
uncompressed n qubits, which belong to half of the n copies of the entangled state \ipg)- Thus, Alice 
and Bob establish n pairs of \ipg). This describes the optimal procedure for the dilution process for 
a pure state. 

The entanglement of distillation and entanglement cost are defined asymptotically, i.e., both 
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processes involve an infinite number of copies of the initial states. For pure states, Eq = Ej) |16j . 
which means that the two processes are reversible asymptotically. Yet, for mixed states, both 
quantities are very difficult to calculate. Nevertheless, it is expected that Ec{p) > Ed{p), viz. that 
one can not distill more entanglement than is put in. 

However, as we now explain, there is a modification of Eq, obtained by averaging Eq over 
pure states, and it is called the entanglement of formation E-p I22j. Any mixed state p can be 
decomposed into mixture of pure states {pi, as in Eq. ()1.6|) . although the decomposition 

is far from unique. To construct the mixed state via mixing pure states in this way will cost, on 
average, ^iPiE{\'ilJi) {tjji\) pairs of Bell states. The entanglement of formation for a mixed state p is 
thus defined as the minimal average number of Bell states needed to realize an ensemble described 
by p, i.e.. 



where the minimization is taken over those probabilities {pi} and pure states that, taken 

together, reproduce the density matrix p = Yl,iPi\'4'i) {'4'i\- Such a construction is usually called a 
convex hull construction. Furthermore, the quantity £'c(|^i)(V'i|) the entropy of entanglement 
of pure state viz. the expression in the right-hand side of Eq. (jl.24p . However, Ep is, in 
general, also difficult to calculate for mixed states, as it involves a minimization over all possible 
decompositions. So far, there has been more analytic progress for E-p than forn Eq and Ep). Notable 
cases include (i) Wootters' formula for arbitary two qubits [22] (or see Appendix[C]); (ii) Terhal and 
Vollbrecht's formula [23j for isotropic states for two qu-dife (d-level parties); and (iii) Vollbrecht 
and Werner's formula [24J for generalized Werner states of two qu-dits. 

One of the central issues in entanglement theory is the so-called additivity of entanglement, 
i.e., whether the entanglement of formation, defined as an average quantity, equals the entangle- 
ment cost, which is defined asymptotically. Recently, Shor [25] has established that the additivity 
problem of entanglement of formation is equivalent to three other additivity problems: the strong 
super additivity of the entanglement of formation, the additivity of the minimum output entropy 
of a quantum channel, and the additivity of the Holevo classical capacity of a quantum channel. 
However, further discussion of these additivity problems is beyond the scope of this dissertation. 




(1.30) 
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1.4 Entanglement via a distance measure 

As any mixture of separable density matrices is still, by definition, separable, any separable state 
can be expressed as a sum of two separable states 

Ps=ppI + {1-p)pI, (1.31) 

unless it is the extremal point, viz. a pure product state. Thus we see that the set of sepa- 
rable states is a convex set. This leads to another type of entanglement measure: the shortest 
"distance" E{p) from an entangled state to the convex set of separable mixed states [26], i.e., 
E{p) = min^eDs d{p\\(T). One example of such an entanglement measure is the relative entropy of 
entanglement, 

£;r(p) EE min Tr(plog/>- plogfj) , (1.32) 
where the distance measure d is defined to be the relative entropy of two states: 

d(p||a) =Tr[plogp-ploga]. (1.33) 

We remark that the relative entropy is non-negative, but it is also not symmetric, i.e., (i(p||(T) 7^ 
d{a\\p). For pure states this definition of entanglement reduces to the entropy of entanglement. 
Another example is the Bures metric of entanglement E^ip), defined via 

Eb{p) = mm[2-2F{p,a)], (1.34) 

where F{p,a) = (Try^ ^fap^faf' is called the fidelity and is symmetric. For two pure states p = 
and cr = the distance d{p\\a) = (2 - 2F{p,a)) reduces to 2(1 - \{^p\(j))\'^). 

We shall give more discussion of the relative entropy of entanglement later on. 

1.5 A simple model 

Let us illustrate how mixed states arise naturally and study the entanglement properties of the 
mixed states. Consider the Hamiltonian H = —Ja^ • ct^. The eigenstates are | jj), | ||), 1^''*') = 
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1 2 3 4 5 



T [fcB/|J|] 

Figure 1.3: Entanglement vs. temperature. The dashed line represents the threshold of entangle- 
ment above which there is a violation of Bell's inequality. 

(I Ti) + I iT))/\/2, and |^'") = (| H) - | iT))/\/2 with respective eigenvalues -J, -J, -J and +3J. 
At a temperature T = 1/(/cb/3) the two-spin system is described by a density matrix 

„I3J p-WJ 1 — r 

P = + + + = + -^t, (1.35) 

where Z = 3e^"^ -t-e"^^*^ and r = (e"^'^'^ — e^'^)/(3e^'^ -j-e"^^"^). This happens to be a Werner state. 
Using the Peres-Horodecki separability criterion, we find that if J > the state is separable, as 
r < 1/3. On the other hand, if J < 0, there is a transition from being entangled to separable as the 
temparature is increased |;29J. This can be seen from the entanglement of formation for the state 
as a function of temperature; see Fig. 11.31 In fact, Werner [llj found that such a state violates the 
Bell-CHSH inequality if r > l/^/2. Thus, there is a finite region between being entangled (r > 1/3) 
and violating the Bell-CHSH inequality. Moreover, Werner found that for r < 1/2 the state can 
be described by a classical local theory. The interval has recently been extended to r < 2/3 by 
Terhal and co-workers [30] . Even more interestingly, as we have seen previously, two pairs of Werner 
states with r > 1/2 can be distilled into a single Werner state of higher r, viz. higher entanglement. 
Thus, given a sufficient supply of Werner states, each of which does not violate any Bell inequality, 
a highly entangled Werner state can be distilled to violate a Bell inequality. This nonlocal quantum 
feature hidden in the initial states is thus revealed by the distillation process [3T] . 

This toy system demonstrates an interesting behavior of entanglement of canonical ensembles 
as the temperature is varied. Later in this dissertation we shall study a more realistic model of 
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large number of spins but at zero temperature. The model that we shall discuss exhibits quantum 
phase transitions as some system parameters vary. 

1.6 Overview of the dissertation 

One of the goals of this thesis is to develop an entanglement theory applicable to many-body (multi- 
partite) systems. But before we apply to more realistic models of many bodies, in Chapter [2] we 
examine the theory in the context of more familiar bi-partite systems, where corresponding results 
for other entanglement measures are already known. We then apply our entanglement measure to 
several nontrivial families of multi-partite states. Connections of our entanglement measure to other 
topics, such as the Hartree approximation, entanglement witnesses, and correlation functions are 
also discussed along the way. In Chapter[3]we discuss in detail the connection between the geometric 
measure of entanglement to the relative entropy of entanglement. In Chapter H] we compute the 
multi-partite measure for two distinct peculiar states, i.e., the so-called bound entangled states. 
Finally, in Chapter [5] we apply our measure to the ground state of the XY quantum spin chain 
model in a transverse magnetic field. The behavior of the entanglement near the quantum critical 
points is found to be dictated by the universality classes of the model. 
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Chapter 2 



Geometric measure of entanglement 
for multi-partite states 

2.1 Introduction 

Only recently, after more than half a century of existence, has the notion of entanglement become 
recognized as central to quantum information processing. As a result, the task of characterizing and 
quantifying entanglement have emerged as one of the prominent themes of quantum information 
theory. There have been many achievements in this direction, primarily in the setting of bi-partite 
systems. Among these, one highlight is Wootters' formula [22] for the entanglement of formation 
for arbitrary two-qubit mixed states. This formula enables discussions of entanglement between 
any pair of two-level systems, which are quite common in various physical systems (or idealizations 
of them). Other achivements include corresponding results for highly symmetrical states of higher- 
dimensional systems \23\ |2l] . 

The success of bi-partite entanglement theories, such as the entanglement of distillation and 
formation, hinges on the reversible interconvertibility of pure entangled states. To be more precise, 
any bi-partite pure entangled states can be, via local operations and classical communication, 
transformed into Bell states, asymptotically and reversibly. This means that there is only one type 
of pure entangled state and one can use Bell states as a "standard ruler" to quantify the degree of 
entanglement. As a result, the entanglement cost Eq equals the entanglement of distillation Ed 
for pure states. 

However, things become more delicate in the multi-partite settings. For example, in the case of 
three qubits, two of the qubits can be entangled whilst the third one is separable from (not entangled 
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with) with them. This kind of entanglement is nothing more than the bi-partite entanglement we 
already know. But there are two other extreme types of entangled states. One is the so-called GHZ 
state: |GHZ) = (|000) + |lll))/^/2, and the other is the W state: \W) = (|001) + |010) + |100))/\/3. 
It was found by Diir and co-workers }32] that, given a single copy of GHZ state, there is no way 
via local operations and classical communication (LOCC), not even probabilistically, that it can be 
transformed into a W state, or vice versa. GHZ and W states have different types of entanglement. 
The GHZ state can be rewritten as 

|GHZ) = i 100)0(1+) + !-)) + ^ |ll)®(|+)-|-)) = l(|OO) + |ll))®|+) + ^(|OO)-|ll))0|-) , (2.1) 

where |ib) = (|0) it \V))/\/2. If one of the parties, e.g., the third one, performs measurement in the 
{lb} basis, and if the outcome is |+), the other two parties have the entangled state l^"*"). If the 
outcome is |— ), the other two parties have the entangled state Therefore, two of the three 

parties can establish a Bell state. Thus, given two copies of GHZ states, they can establish a Bell 
state shared between say A and B, and another Bell state shared between A and C. A can locally 
prepares a W state 

|W) = i=(|0i02l3) + IO1I2O3) + IliOaOg)). (2.2) 

As she has one Bell state shared with B and the other with C, she can use quantum teleportation 
to teleport the state of particle 2 of the W state to B and that of particle 3 to C. Thus two copies 
of GHZ states can achieve a copy of the W state. 

On the other hand, if the three parties share a W state and one of them makes a measurement 
in the {0/1} basis, 2/3 of the time the other two parties can establish a Bell state. The other 1/3, 
they fail to do so. If they are given 2 copies of W states, 4/9 of the time they can establish one Bell 
state between, say, A and B, and the other Bell state between A and C. A can use the same trick 
of teleportation such that they end up with a GHZ state. So two copies of W states can achieve 
4/9 copy of GHZ state. 

However, it is still not yet clear whether or not given some number of copies of GHZ states, 
they can be transformed into as many copies of W states and then transformed back to the original 
number of GHZ states. Namely, it is not yet known that whether the process is asymptotically 
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reversible. 

Bennett and co-workers have come up with the notion of a finite minimal reversible entangle- 
ment generating set (MREGS) |33j- Such a set should include different types of "ruler" states. 
But the numer of states in this set should be finite, otherwise, it is not practical to use an in- 
finite number of different rulers to "measure" the property of entanglement. For example, the 
minimal set might contain { | Bell) -|^2 1 1 I^^^l) 13 ' 1^^611)23 , | GHZ)}. Perhaps a large number of any 
three-qubit pure state can be transformed reversibly into, say, xi copies of |Bell)j^2i ^2 copies of 
I Bell) -1^3, X3 copies of | Bell) 23, and X4 copies of |GHZ)? Or perhaps the minimal set might contain 
{|Bell)-^2 ; 1^^^11)13 5 11^611)23 ' |GHZ) , |VF)}? Then the entanglement could be defined as some kind 
of vector. However, this MREGS problem has not been solved yet. The situation gets even worse 
beyond three qubits. Verstraete and co-workers [Mj have found that there are nine inequivalent 
classes of four-qubit entangled states. As the dimensions and the number of parties grow, the 
number of states in MREGS, if such a set exists, is expected to grow considerably large. 

The above considerations complicate the task of extending measures such as the entanglement 
of distillation [16j and formation |201 l22j to multi-partite systems. Moreover, we have seen that the 
characterization of general multi-partite entanglement remains incomplete, as the number of the 
types of entanglement grows with the number of parties and the dimensions of Hilbert space. The 
issue of entanglement for multi-partite states hence poses an even greater challenge than bi-partite 
states. On the other hand, one can quantify multi-partite entanglement via other measures, such 
as the relative entropy of entanglement, the Bures metric [26^ I27j. and the Schmidt measure |35j . 
which are naturally extendible to multi-partite settings. 

In this chapter, we present an attempt to quantify multi-partite entanglement by developing and 
investigating a certain geometric measure of entanglement (GME), first introduced by Shimony [36] 
in the setting of bi-partite pure states and generalized to the multi-partite setting (via projection 
operators of various ranks) by Barnum and Linden [37J. In Sec. 12.21 we begin by examining this 
geometric measure in pure-state settings and establishing a connection with entanglement witnesses, 
Hartree approximations, and correlation functions. In Sec. 12.31 we extend the measure to mixed 
states, showing that it satisfies certain criteria required of good entanglement measures. In Sec. 12.41 
we examine the GME for several families of mixed states of bi-partite systems: (i) arbitrary two- 
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Figure 2.1: The schematic picture of the geometric measure. Imagine all pure states lie on the 
sphere and all separable pure states |(/))'s lie on the equator. The degree of entanglement for is 
reflected in the shortest distance to the set of separable states. 

qubit mixed, (ii) generalized Werner, and (iii) isotropic states in bi-partite systems, as well as 
(iv) certain mixtures of multi-partite symmetric states. In Sec. 12.51 we give a detailed application 
of the GME to arbitrary mixtures of three-qubit GHZ, W and inverted-W states. In Sec. 12.61 we 
discuss some open questions and further directions. The discussion in this chapter is based on 
Ref. [3H] 

It is not our intention to cast aspersions on existing approaches to entanglement; rather we 
simply wish to add one further element to the discussion. Our discussion focuses on quantifying 
multi-partite entanglement rather than characterizing it. 

2.2 Basic geometric ideas and application to pure states 

We begin with an examination of entangled pure states, and of how one might quantify their 
entanglement by making use of simple ideas of Hilbert space geometry. Let us start by developing 
a quite general formulation, appropriate for multi-partite systems comprising n parts, in which 
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each part can have a distinct Hilbert space. Consider a general n-partite pure state 



1^^)= E w-pJeff4?---4"3' (2-3) 

Pl---Pn 

where {cp^^} is the local basis of the k-th party, e.g., the spin | or |. One can envisage a geometric 
definition of its entanglement content via the distance 

d = min|| 1^) - |(/))|| (2.4) 
10} 

between {tp) and the nearest of the separable states {(p) (or equivalently the angle between them). 
Here 

= = \^^^^) • • • (2.5) 

is an arbitrary separable (i.e., Hartree) n-partite pure state, the index i = 1 . . .n labels the parties, 
and a state vector of part i is written as 



It seems natural to assert that the more entangled a state is, the further away it will be from its best 
unentangled approximant (and, correspondingly, the wider will be the angle between them). We 
emphasize that we only compare the entangled pure state to the set of pure unentangled state. We 
shall extend to mixed states via the so-called convex-hull construction in Sec. l2.3i Another approach 
one might take is to compare any entangled state to the set of unentalged states, including both 
pure and mixed. The Bures measure, introduced in Sec. II. 4^ is such an example. 

To actually find the nearest separable state, it is convenient to minimize, instead of d, the 
quantity d?, i.e., 

(2.7) 



subject to the constraint = 1. In fact, in solving the resulting stationarity condition one 

may restrict one's attention to the subset of solutions \<j)) that obey the further condition that each 
factor obeys its own normalization condition = 1. Thus, by introducing a Lagrange 
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mulitplier A to enforce the constraint {4>\4') = 1, differentiating with respect to the independent 
ampUtudes, and then imposing the further condition = 1, one arrives at the nonlinear 

eigenproblem for the stationary 

E x;„„c« •••c'« •••cW = A4f , (2.8a) 
V V . . . ^W* . . . - Ar(*) 1^2 8b) 



PVPVPn 



where the eigenvalue A is associated with the Lagrange multiplier enforcing the constraint = 1, 
and the symbol ^ denotes the exclusion of the corresponding term or factor. In a form independent 
of the choice of basis within each party, Eqs. (j2.8p read 



,\( \cp(j))) = a/(A« , (2.9a) 
{cj)^^^)\i;) = A (2.9b) 

From Eqs. (12. 8p or (12. 9p . e.g., by taking inner product of both sides of Eq. (I2.9ap with |(/>^*^) one 
readily sees that 

K={m = {m (2.10) 

and thus the eigenvalues A are real, in [—1,1], and independent of the choice of the local basis 
(kp?)}- Hence, the spectrum A is the cosine of the angle between and \4>)\ the largest, Amax, 
which we call the entanglement eigenvalue, corresponds to the closest separable state and is equal 
to the maximal overlap 

A^ax = max||(0|^)||, (2.11) 
<t> 

where \(j)) is an arbitrary separable pure state. 

Although, in determining the closest separable state, we have used the squared distance between 
the states, there are alternative (basis- independent) candidates for entanglement measures which 
are related to it in an elementary way: the distance, the sine, or the sine squared of the angle 
9 between them (with cos6' = Re (-010)). We shall adopt -Egjjj2 = 1 — A^^^^ as our entanglement 
measure because, as we shall see, when generalizing E^:^^2 to mixed states we have been able to 
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show that it satisfies a set of criteria demanded of entanglement measures. 

In bi-partite apphcations, the eigenproblem (12. Sp is in fact linear, and solving it is actually 
equivalent to finding the Schmidt decomposition [36j. Moreover, the entanglement eigenvalue is 
equal to the maximal Schmidt coefficient. To be more precise, in bi-partite systems the stationarity 
conditions (j2.8p reduce to the linear form 

E4iP.4'? =A4f, (2.12a) 
Pi 

Pi 

E^;iP.4? =^4f' (2-12C) 

P2 

^XpiP24f =^4i^- (2-12d) 



P2 



Eliminating c^^ between Eqs. (j2.12ap and (j2.12dp and, similarly, eliminating c^^ between Eqs. (|2.12bp 
and (j2.12cp gives 

T.^P^P^K[P2''i! = ^'^ff' (2.13a) 



PiP2 



^PiP2 Pa 



P1P2 



or equivalently 



Tn(\^MW'^) = A2|0(2)). (2.14b) 



Now, solving the above equations is equivalent to finding the Schmidt decomposition for {ip). 
To see this, first recall that Schmidt decomposability guarantees that an arbitrary pure bi-partite 
state 

|V') = 2^Xp,p2|4i^)»|420 (2-15) 
P1P2 

in the (product) Hilbert space 7ii®7i2 (with factor dimensions di and ^2) can always be expressed 
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in the form 

min{di,d2} 



E ^> 4^0- (2-16) 



2(2) 



k=l 

Here, Afc > 0, J2k^k ~ ^'^^ I^^V^ and |e[,^'*)'s are orthonormal, respectively, in Tii and 7^2- 
Moreover, the new (tilde) bases, as well as the A^'s, follow as the solution of the eigenproblems of 
the reduced density matrix that one obtains by tracing [tp] over party 1 or 2: 

^2(1^X^1)14'^) = A^l4'^)' (2-17a) 
Tri(|V')(^l)l4'^) = Ailef)). (2.17b) 

These are Eqs. (j2.14p : thus we see that determining entanglement for bi-partite pure states is 
equivalent to finding their Schmidt decomposition, except that one only needs the largest Schmidt 

coefficient Amax = Amax- 

By contrast, for the case of three or more parts, the eigenproblem is a nonlinear one. For 
example, in the setting of tri-partite systems, the stationarity conditions (|2.8|) associated with the 
pure state l^-) = EpipaPs ^piP2P3|41^42''43^) become 



/ ^ ^PlP2P3 P2 ^P'i 




P2P3 




/ J -^PlP2P3 Pl P3 


= Ac(2)* 


PlP3 




/ J ■^PlP2P3 Pi P2 


= Ac^^)* 


P1P2 





(2.18a) 
(2.18b) 
(2.18c) 



Note the nonlinear structure of this tri-partite (and, in general, any n > 3-partite) eigenproblem. 
As far as we are aware, for nonlinear eigenproblems such as these, one cannot take advantage 
of the simplicity of linear eigenproblems, for which one can address the task of determining the 
eigenvalues directly (via the characteristic polynomial), without having to address the eigenvectors. 
Hence, even for systems comprising qubits, one is forced to proceed numerically. This is consistent 
with the notion that no Schmidt decomposition exists beyond bi-partite systems. Yet, as we shall 
illustrate shortly, there do exist certain families of pure states whose entanglement eigenvalues can 
be determined analytically. 
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For C"^ (8) bi-partite systems, the equivalence between the geometric approach and Schmidt 
decomposition immediately indicates why the maximally entangled pure states have (up to local 
unitary transformations) the well-known form \ ^^) = Yli=i N)- Schmidt coefficients are 

non- negative and sum (when squared) to unity, any less symmetric state must have a larger Amax; 
i.e., a smaller entanglement. 

2.2.1 Illustrative examples 

Suppose we are already in possession of the Schmidt decomposition of some two-qubit pure state: 



1^) = VP|00) + Vl^lll)- (2-19) 
Then we can read off the entanglement eigenvalue: 

Amax = max{^, a/1 -p}. (2.20) 



Another approach to obtain A^ax is to solve, e.g., Eqs. (|2.14p . which lead to solving the maximal 
eigenvalue of 

/ 

p 



v = A^v, (2.21) 



1-p^ 



resulting in the solution (j2.20p . 



Now, recall [22] that the concurrence C for this state is 2y/p(l — p) [see Appendix [C]. Hence, 
one has 

ALx = ^ (l + , (2.22) 

a connection between the entanglement eigenvalue and the concurrence that holds for arbitrary 
two-qubit pure states. 

The possession of symmetry by a state can alleviate the difficulty associated with solving the 
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nonlinear eigenproblem. To see this, consider a state 



'lP2--Pn 




(2.23) 



Pl—Pn 



that possesses the symmetry that the (nonzero) amphtudes x ^-re invariant under permutations. 
What we mean by this is that, regardless of the dimensions of the factor Hilbert spaces, the am- 
plitudes are only nonzero when the indices take on the first z/ values (or can be arranged to do 
so by appropriate relabeling of the basis in each factor) and, moreover, that these amplitudes are 
invariant under permutations of the parties, i.e., Xo-io-2--o-„ = XpiP2-Pni where the a's are any per- 
mutation of the p's. (This symmetry may be obscured by arbitrary local unitary transformations.) 
For such states, it seems reasonable to anticipate that the closest Hartree approximant retains this 
permutation symmetry. Assuming this to be the case — and numerical experiments of ours support 
this assumption — in the task of determining the entanglement eigenvalue one can start with the 
Ansatz that the closest separable state has the form 



i.e., is expressed in terms of copies of a single factor state, for which Cj = Cj. To obtain the 
entanglement eigenvalue it is thus only necessary to maximize Re with respect to the v 

amplitudes {cj}^^^, a simpler task than maximization over the X^iLi amplitudes of a generic 
product state. 

To illustrate this symmetry-induced simplification, we consider several examples involving 
permutation-invariant states, first restricting our attention to the case v = 2 (i.e., two-state par- 
ties). The most natural realizations are n-qubit systems. One can label these symmetric states 
according to the number of O's, as follows [39j : 



permutations k n—k 

As the amplitudes in this state are all positive and the state is permutationally invariant, one can 
assume that the closest Hartree state also has the symmetry and has non-negative amplitudes, and 




(2.24) 




(2.25) 
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is of the form 

!</-) = (V^|0) + yr^ 1 1))^" (2.26) 
Maximizing the overlap with respect to p, one obtains the entanglement eigenvalue for \S{n,k)): 



k_ n—k 

k\ 2 / n - k\ 2 



For fixed n, the minimum Amax (and hence the maximum entanglement) among the |S(n, fe))'s 
occurs for k = n/2 (for n even) and k = {n i: l)/2 (for n odd). In fact, for fixed n the general 
permutation-invariant state can be expressed as J2k ^)) "'^ith Y^j^ la^p = 1. The entangle- 

ment of such states can be addressed via the strategy that we have been discussing, i.e., via the 
maximization of a function of (at most) three real parameters. The simplest example is provided 
by the nGHZ state: 

InGHZ) = ^ ( |S(n, 0)) + |S(n, n)) ) . (2.28) 
v2 

It is easy to show that (for all n) Ainax(nGHZ) = l/^/2 and, equivalently, E^^^2 = 1/2. Note that 
one could have rescale the definion of the GME by a factor of two such that Bell and N-GHZ states 
have entanglement of unity. However, we do not make this rescaling. 

We now focus our attention on three-qubit settings. Of these, the states |5(3,0)) = |000) 
and |5(3,3)) = |111) arc not entangled and are, respectively, the components of the 3-GIIZ state: 
|GHZ) = ( |000) + |111))/V2. The states |S(3,2)) and |S(3, 1)), are the W state and the "inverted 
W" state, respectively, 

|W) = |S(3,2)) = -^(|001) + |010) + |100)), (2.29a) 
W) ^ |S(3,1)) = -^(|110) + |101) + |011)), (2.29b) 

are equally entangled, having Amax = 2/3 and ^ginZ = 1 — A^g^ = 5/9. 
Next, consider a superposition of the W and W states: 



WW(s,(/))) = ^|W) + ^r^e''^ W). (2.30) 
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Figure 2.2: Entanglement of the pure state ^/s\W) + \/l — s |W) vs. s. This also turns out to be 
the entanglement curve for the mixed state s |W)(W| + (1 — s)|W)(W|. 



It is straightforward to see that its entanglement is independent of (/>: the transformation | |0) , |1) } — 
{ |0) jB"**^ |1) } induces WW(s,(/>)^ e"**^ WW(s,0)y and the entanglement is independent of 
the overall phase of a state. To calculate Amax) we assume that the separable state is of the form 
(cos 9 |0) + s'mO and maximize its overlap with WW(s, 0)^ with respect to 9. Thus we find 

that 9{s) is determined via its tangent t = tan^, which is the particular root of the polynomial 
equation 

Vl-st^ + 2^5 - 2Vl - s t - ^/s = (2.31) 

that hes in the range t G [y^l/2, a/2]. Via 9{s), Amax (and thus E^-^j^2 = 1 — A^^^) can be expressed 
as 



Amax(s) = -^[a/scos0(s)+ \/ 1 — s siu 9 {s)] sin26'(s). 



(2.32) 



In Fig. [221 we show £^35^2 ( WW(s, ) vs. s. 
In fact, for the more general superposition 



= \Sin,h)) + Vl^e"^ |S(n,A:2)) 



(2.33) 



(with ki ^ k2) A 

max also turns out to be independent of (p, and can be computed in the same 



way as was used for 



WW(s, (p) ). We note that although the curve in Fig. 12.21 is convex, convexity 
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Figure 2.3: Entanglement of I^ghz+wI-Sj 4>)) vs. s. The upper curve is for (/> = tt whereas the lower 
one is for (/> = 0. Dots represent states with randomly generated s and (p. 

does not hold uniformly over ki and k2- This adds some complexity to calculation of mixed-state 
entanglement, as we shall see later. 

As our last pure-state example in the qubit setting, we consider superpositions of W and GHZ 
states: 



For these, the phase (p cannot be "gauged" away and, hence, E^:^^2 depends on (j). In Fig. 12.31 we 
show -Bsin^ vs. s at </) = and cj) = tt (i.e., the bounding curves), as well as E^^^2 for randomly 
generated values of s G [0,1] and (p G [0, 27r] (dots). We observe that the 'vr' state has higher 
entanglement than the '0' does. As the numerical results suggest, the (</>-parametrized) E^^^2 vs. s 
curves of the states |^ghz+w('5, (/>)) lie between the 'vr' and '0' curves. 

We remark that, more generally, for systems comprising n parts, each being a d-level system, 
the symmetric state 



|^ghz+w( 



s, (f))) = ^fs I GHZ) + y/l^s e'^ |W) . 



(2.34) 




(2.35) 



(with ki = n) has entanglement eigenvalue 




(2.36) 
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One can also consider other symmetries. Consider the totally anii-symmetric (viz. determinant) 
state of n parts, each with n levels, 

1 " 

|Det„) EE — = V eii,,.,,i„ . . . • (2.37) 

Vn! . "—^ , 

Ji,...,i„=i 

It has been shown by Bravyi [30] that the maximal squared overlap of this state A^^^^ is 1/nl. 
Bravyi has also generalized the anti-symmetric state to the n = p (i^-partite determinant state, via 
the construction 

(/.(I) = (0,0,..., 0,0), 
(l){2) = (0,0,..., 0,1), 

(j){dP-l) = {d-l,d-l,... ,d-l,d-2), 
^{dP) = {d-l,d-l,...,d-l,d-l), 

and 

|Det„,d) = J— J2 ^ii,-,idp\(l^ih),---,<P{idp)) ■ (2.38) 
In this case, we can show that A^^^^ = [(^^''')']^^- 
2.2.2 Connection with the Hartree approximation 

Recall from Eq. (j2.1ip that the maximal overlap or the entanglement eigenvalue of a pure state 
1^) is defined as Ajnax(V') = max,^ where the maximization is taken over arbitrary product 

state 10). On the other hand, consider the Hamiltonian constructed from the entangled state {ip) 

H^ = -|V)(VI, (2.39) 

which has the form of a projector on to the state {ip) and which has the minimum eigen-energy — 1. 
The Hartree approximation is to take some product state \^) such that it minimzes the expectation 
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value of the Hamiltonian 



min(0| W^|</.) = -max|(,/.|V')P = -ALxW- (2-40) 

9 <P 

It thus arises that the difference between the true ground-state energy of TC^ and the one given by 
the best Hartree approximation (j2.40p is 

1-ALxW, (2.41) 

i.e., the GME of {ip) and that that the best Hartree approximant {(p*) is the closest separable state 
to IV^). 

This connection, although elementary, can be generalized to the Hartree-Fock approximation 
when we wish to extend the idea of entanglement to fermionic systems. 

If we are interested in the entanglement of a nondegenerate ground state (with energy £"0) of a 
Hamiltonian, in principle, we can express the maximal overlap as 

ALx(|i^o» = lim min (<A| e-'^^^"^") |0) . (2.42) 

2.2.3 Connection with entanglement witnesses 

We now digress to discuss a relationship between the geometric measure of entanglement and 
another entanglement property — entanglement witnesses. The entanglement witness W for an 
entangled state p is defined to be an operator that is (a) Hermitian and (b) obeys the following 
conditions [13]: 

(i) Tr(Wcr) > for all separable states cr, and 

(ii) Tr(>Vp) < 0. 

Here, we wish to establish a correspondence between Amax for the entangled pure state IV') and the 
optimal element of the set of entanglement witnesses W for that have the specific form 

W = A^l - (2.43) 
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where 1 is the identity operator. This set is parametrized by the real, non-negative number A^. By 
optimal we mean that, for this specific form of witnesses, the value of the expression Tr(W|V')(V'|) 
is as negative as can be. 

In order to satisfy condition (i) we must ensure that, for any separable state a, we have 
Tr(W(T) > 0. As the density matrix for any separable state can be decomposed into a mixture 
of separable pure states [i.e., a = J2i where {|</>i)} are separable pure states], condition (i) 

will be satisfied as long as Tr(>V|(/)) (01) > for all separable pure states This condition is 
equivalent to 

- IKV'I'^)!!^ > (for all separable | (2.44) 

which leads to 

\''>maMm' = ^L.m)- (2.45) 

m 

Condition (ii) requires that Tr(yV|^)(^|) < 0, in order for W to be a valid entanglement witness 
for 1^); this gives - 1 < 0. Thus, we have established the range of A for which A^l — is a 

valid entanglement witness for 

ALx(IV')) < A2 < 1. (2.46) 

With these preliminaries in place, we can now establish the connection we are seeking. Of the 
specific family (I2.43P of entanglement witnesses for the one of the form Wopt = ^maxdV'))! ~ 
\ip){ip\ is optimal, in the sense that it achieves the most negative value for the detector Tr (Woptl V') (^1) ^ 

niinTr(W|V)(V'l) = IVC^Voptl^X^I) = -i^sin^d^)), (2.47) 

where W runs over the class (j2.43p of witnesses. 

We now look at some examples. For the GHZ state the optimal witness is 

Wghz = ^1-|GHZ)(GHZ| (2.48) 

and Tr(WGHz|GHZ)(GHZ|) = -£^35^2 (| GHZ)) = -1/2. Similarly, for the W and inverted-W states 



32 



we have 

Ww = ^l-|W)(W| and >V^ = ^1-|W)(W| (2.49) 
and Tr(>Vw|W)(W|) = -E^^^2{\W)) = -5/9, and similarly for For the four-qubit state 

l^-) = (|0011) + |0101) + |0110) + |1001) + |1010) + |1100))/a/6 (2.50) 

the optimal witness is 

3 

Wvi, = -1 - |^')(^| (2.51) 
8 

and Tr(W^|^)(^'|) = -S,i„2(|*)) = -5/8. 

Although the observations we have made in this section are, from a technical standpoint, elemen- 
tary, we nevertheless find it intriguing that two distinct aspects of entanglement — the geometric 
measure of entanglement and entanglement witnesses — are so closely related. Furthermore, this 
connection sheds new light on the content of the geometric measure of entanglement. In particular, 
as entanglement witnesses are Hermitian operators, they can, at least in principle, be realized and 
measured locally }41j . Their connection with the geometric measure of entanglement ensures that 
the geometric measure of entanglement can, at least in principle, be verified experimentally. 

2.2.4 Connection with correlation functions 

In this section we explore the connection of the entanglement eigenvalue to the correlation functions. 
To illustrate the connection, we shall mainly focus on A^-qubit systems. Recall that a single pure 
qubit state, when written in the form of a density matrix, can be expressed in terms of Pauli 
matrices plus the identity: 

^(l + f-a), (2.52) 

with the constraint on the real vector parameter r that |r| = 1. The density matrix for an A^-qubit 
separable pure state {(p) ()2.5p is hence a direct product of A^ such terms: 

\4>){4>\ = i(l + fi • ffW) ^(1 + f2 ■ a^^)) • • • ® ^(1 + rjv • <?^^^), (2.53) 
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and expanding out the products gives the decomposition of \4'){4'\ 




N 



j=l j<k 



j<k<l 



<J) ffc • (t(*^) f; • a(') + • • • }. (2.54) 



The overlap squared of any A'"-qubit state with the separable state \(p) is 



A2 



(2.55) 



which can in turn be expressed in terms of the correlations of the operators of the form a (8> • • • (8) 
a <Si ■ ■ ■ ■ Therefore, the maximal overlap is a maximization over all combinations of M-point (with 
M < N) correlation functions: 



The average is taken with respect to the A^-qubit pure state I?/') whose entanglement we are to 
quantify. This result echos the fact mentioned in the previous section that the geometric measure 
(or here the maximal overlap) can be measured, via the correlations, which are physical quantities. 

2.3 Extension to mixed states 

The extension of the geometric measure to mixed states p can be made via the use of the convex 
roof (or hull) construction [indicated by "co"], as was done for the entanglement of formation [c.f. 
Eq. (|1.30p ] [22] . The essence is a minimization over all decompositions p = YliPi into pure 

states, i.e.. 



This convex hull construction ensures that the measure gives zero for separable states; however, in 
general it also complicates the task of determining mixed-state entanglement. 




(2.56) 




(2.57) 
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Now, any good entanglement measure E should, at least, satisfy the following criteria (c.f. Refs. \26\ 

Cl. (a) E{p)>0; (b) E{p) = if /? is not entangled. 
C2. Local unitary transformations do not change E. 

C3. Local operations and classical communication (LOCC) (as well as post-selection) do not 
increase the expectation value of E 0. 

C4. Entanglement is convex under the discarding of information, i.e., Y2iPi ^(Pi) — ^i^iPi Pi)- 

The first requirement simply states that the entanglement is a non-negative quantity, and it vanishes 
for unentangled states. Local unitary transformations are simply a change of local basis, and 
entanglement should be invariant under such a change. The third criterion simply requires that 
the average entanglement cannot be increased during manipulations that are local, which reflects 
the fact that entanglement is a nonlocal resource. The last criterion states the fact that, for a set 
of states, entanglement can never be increased if the information about which state is which is 
lost. The issue of the desirability of additional features, such as continuity and additivity, requires 
further investigation, but C1-C4 are regarded as the minimal set. If a measure satisfies C1-C4, it 
is called an entanglement monotone [42]. 

Does the geometric measure of entanglement obey Cl-4? The answer is affirmative but we shall 
relegate the proof to Appendix |Dj We remark that from the definition (j2.57p it is evident that Cl 
and C2 are satisfied, provided that -Epure satisfies them. Furthermore, the convex-hull construction 
automatically fulfills C4. The consideration of C3 seems to be more delicate and whether or not 
it holds depends on the explicit form of impure- In Appendix |D] we show that the choice of taking 
£^gj^2 as the entanglement measure does satisfy C3. 

2.4 Analytic results for mixed states 

Before moving on to the terra incognita of mixed multi-partite entanglement, we test the geometric 
approach in the setting of mixed bi-partite states, by computing E^^^2 for three classes of states for 
which Ep is known. 

^This requirement does not contradict with distillation, as it takes into account the cases when distillation fails. 
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2.4.1 Arbitrary two-qubit mixed states 



For arbitrary two-qubit states we can show that 



(2.58) 



where C{p) is the Wootters concurrence of the state p. The essential point is that, in his derivation of 
the formula for E-p^ Wootters showed that there exists an optimal decomposition p = ^^Pi \ilJi){ipi\ 
in which every {ipi) has the concurrence of p itself. (More explicitly, every has the identical 
concurrence, that concurrence being the infimum over all decompositions.) By using Eq. (I2.22[) we 
see that Eq. ()2.58p holds for any two-qubit pure state. As -Esin^ is a monotonically increasing convex 
function of C G [0, 1], the optimal decomposition for E^^^2 is identical to that for the entanglement 
of formation Ep. Thus, we see that Eq. (j2.58p holds for any two-qubit mixed state. 

The fact that E^^^2 is related to Ep via the concurrence C is inevitable for two-qubit systems, 
as both are fully determined by the one independent Schmidt coefficient. We note that Vidal [44J 
has derived this expression when he considered the probability of success for converting a single 
copy of some pure state into the desired mixed state p. The probability of the conversion is 



which gives a physical interpretation of the geometric measure of entanglement. Unfortunately, 
this connection only holds for two-qubit states. 

2.4.2 Generalized Werner states 

Any state of a (E> C"^ system is called a generalized Werner state if it is invariant under 



where U is any element of the unitary group U{d) and dU is the corresponding normalized Haar 
measure. Both parties perform the same unitary transformation U on their own system but they 
choose U randomly according to the measure dU. 




(2.59) 




(2.60) 
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Generalized Werner states can be expressed as a linear combination of two operators: the 
identity 1, and the swap F = - \ij){ji\, i-e., pw = al + ^F, where a and b are real parameters 
related via the constraint Trpw = 1- This one-parameter family of states can be neatly expressed 
in terms of the single parameter / = Tr(/9'vvF): 

^,^. + Z^f. ,2.61) 

By applying to E^^^2 the technique by developed by Vollbrecht and Werner for i5^F(pw) [to be briefly 
reviewed shortly; or see Ref. [24j], one arrives at the geometric entanglement function for Werner 
states: 

^sin2(pw(/)) = ^(l-\/l-/2) for/<0, (2.62) 

and zero otherwise. 

In order to prepare for this, we now briefly review the technique developed by Vollbrecht and 
Werner [23] for computing the entanglement of formationfor pw ; this turns out to be applicable to 
the computation of the sought quantity E^^^2 . We start by fixing some notation. Let 
i^' be a compact convex set (e.g., a set of states that includes both pure and mixed ones); 
M be a convex subset of K (e.g., set of pure states); 

E : M ^ RU {+00} be function that maps elements of M to the real numbers (e.g., E = E^^^2); 
and 

G be a compact group of symmetries, acting on K (e.g., the group U x W) as ag : K ^ K (where 
Ug is the representation of the element g ^ G) that preserve convex combinations. 
We assume that OgM C M (e.g., pures states are mapped into pure states), and that E{agm) = 
E{m) for all m G M and g £ G (e.g., that the entanglement of a pure state is preserved under ag). 
We denote by P the invariant projection operator defined via 

P/c = y dgag{k), (2.63) 

where k £ K. One example of Pi is the operation in Eq. (j2.60p . Vollbrecht and Werner also defined 
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the following real-valued function e on the invariant subset PK: 



e{x) = inf {E{m)\m G M, Pm = x} . 



(2.64) 



They then showed that, for x G PK, 



CO E(x) 



CO e{x) 



(2.65) 



and provided the following recipe for computing the function co E for G-invariant states: 

1. For every invarinat state p (i.e., obeying p = Pp), find the set Mp of pure states a such that 
Pa = p. 

2. Compute e(p) = inf {£;(cr)|cr € Mp}. 

3. Then co E is the convex hull of this function e. 

Having reviewed the Vollbrecht- Werner technique, we now apply it to the geometric measure 
of entanglement E^^^2 by applying their recipe. 

The essential points of the derivation are as follows: 
(i) In order to find the set Mp^ it is sufficient, due to the invariancc of pw under Pi, to consider 
any pure state |$) = J2jk^jk\^f^) ki^'*) ^^^^ ^ diagonal reduced density matrix Tr2|$)($| 
and the value Tr(|$)($|F) equal to the parameter / associated with the Werner state pw(/)- It 
can be shown that 



where Xu = |$m| . 

(ii) If / > 0, we can set the only nonzero elements of |$) to be ^n, ^i2, ■ ■ ■, ^u, • • •, ^id such that 
l^iiP = /) this state obviously being separable. Hence, for / > we have £^g;jj2(pw(/)) = 0. On 
the other hand, if / < then any nonzero Xu would increase (/ — J2i ^u)^ a^nd, hence, increase the 
value of £?(!$) ($1), not conforming with the convex hull. Thus, for a fixed value of /, the lowest 
possible value of the entanglement E{\^){^\) that can be achieved occurs when Xu = and there 




(2.66) 
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are only two nonzero elements and (i 7^ j). This leads to 

E{\^)m = ^[i-y^r^). (2.67) 

Thus, as a function of /, e(/) is given by 



mm 

|<I>) at fixed / 



i l-TW^ for/<0, 
e(/) = <J ' ^ ^ (2.68) 

for / > 0, 

which, being convex for / € [—1, 1], gives the entanglement function (12.62P for Werner states. 

2.4.3 Isotropic states 

Instead of both performing U on their system, one of the two parties performs U whereas the other 
performs U* . The isotropic states are are invariant under 

P2- J dU{U0U*)p{U^ 0U*^), (2.69) 

and can be expressed as 

Aso(F) = ^2=^ (i - If +)(<J>+|) + i^|f +)(f +1, (2.70) 

where |^>+) = Y^i=i and F G [0, 1]. For F £ [0, 1/d] this state is known to be separable 03]. 
By following arguments similar to those applied by Terhal and Vollbrecht [23j for £'f(Aso) one 
arrives at 

E,,^. (Aso(F)) = 1 - ^(Vf+ V(l-F)(d-l))^ (2.71) 

for F > 1/d. The essential point of the derivation is the following Lemma (c.f. Ref. |23j): 
Lemma 1. The entanglement E^^^2 for isotropic states in (E) for F G [1/d, 1] is given by 

E^,^2MF)) = coiR{F)), (2.72) 



39 



0.68 
. 67 
.66 
. 65 

0.2 0.4 0.6 0.8 1 

r 

Figure 2.4: Entanglement curve for the mixed state P7-2,5i''') (full line) constructed as the convex 
hull of the curve for the pure state |SS7;2,5(?', cp)) (dashed in the middle; full at the edges). 

where co{R{F)) is the convex hull of the function R and 

d 2 ^ 

i?(F) = l-max{/ii|F= (VV^) /d; V/Xi = l}. (2.73) 
Straightforward extremization shows that 



R{F) = 1 - 




(2.74) 



which is convex, and hence co{R{F)) = R{F). Thus we arrive at the entanglement result for 
isotropic states given in Eq. ()2.7ip . 

2.4.4 Mixtures of multi-partite symmetric states 

Before exploring more general mixed states, it is useful to first examine states with high symmetry. 
With this in mind, we consider states formed by mixing two distinct symmetric states (i.e., ki / k2)- 

Pn-Mk,{r) = r |S(n, A:i)) (S(n, k,)\ + (1 - r) |S(n, ^2)) (S(n, A^a)! . (2.75) 
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From the independence of E^^^2 {\SSn;kik2{^T (p))) o^i and the fact that the mixed state Pn-Mkjy^) 
is invariant under the projection 

P3 : [ ^u^n^^i^n (2.76) 

J 27r 

with ?7 : { |0) , |1) } ^ { |0) ,e~"^ |1) }, we have that E^^^2 (Pn;fciA:2 (^)) be constructed from 

the convex hull of the entanglement function of |SS„;fcifc2(^i 0)) example, {n,ki,k2) = 

(7, 2, 5), is shown in Fig. 12.41 If the dependence of -^81^2 on r is already convex for the pure state, 
its mixed-state counterpart has precisely the same dependence. Figure [2T2l for which {n,ki,k2) = 
(3, 1, 2), exemplifies such behavior. More generally, one can consider mixed states of the form 

P(M) = (2.77) 
k 

The entanglement -Emixed({p}) can then be obtained as a function of the mixture {p} from the 
convex hull of the entanglement function Ep^^(.({q}) for the pure state Ylk ^))' That is, 

Emixiip}) = coii^pure({^})|{g=p}- Therefore, the entanglement for a mixture of symmetric states 
\S{n,k)) is known from Ep^j-c{{q}), up to some convexification. 

2.5 Application to arbitrary mixture of GHZ, W and inverted-W 
states 

Having warmed up in Sec. 12.4.41 by analyzing mixtures of multi-partite symmetric states, we now 
turn our attention to mixtures of three-qubit GHZ, W and inverted-W states. 

2.5.1 Symmetry and entanglement preliminaries 

These states — GHZ, W and inverted-W states — are important, in the sense that all pure three-qubit 
entangled states can, under probabilistic LOCC, be transformed either to GHZ or W (equivalently 
inverted-W) states. It is thus interesting to determine the entanglement content (using any measure 
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of entanglement) for mixed states of the form: 



p(x,y) =x|GHZ)(GHZ| + y|W)(W| + (l-a;-y)|W)(W|, (2.78) 

where x,y > and x + y < 1. This family of mixed states is not contained in the family (j2.77p . 
as |GHZ) = (|5(3,0)) + \S{3,3)))/\/^. The property of p{x,y) that facilitates the computation 
of its entanglement is a certain invariance, which we now describe. Consider the local unitary 
transformation on a single qubit: 

|0) ^ |0), (2.79a) 
|1) ^ (2.79b) 

where g = exp (27ri/3), i.e., a relative phase shift. This transformation, when applied simultaneously 
to all three qubits, is denoted by Uk- It is straightforward to see that p{x, y) is invariant under the 
mapping 

1 ^ 

^A-.p^-Y^UkpUl (2.80) 

Thus, we can apply Vollbrecht- Werner technique [23] to the compution of the entanglement of 
p{x,y)- 

Now, the Vollbrecht- Werner procedure requires one to characterize the set Si-av of all pure states 
that are invariant under the projection P4. Then, the convex hull of E^^^2{p) need only be taken 
over Sinv, instead of the set of all pure states. However, as the state p{x,y) is a mixture of three 
orthogonal pure states (viz. |GHZ), |W) and W^) that are themselves invariant under P4, the 
pure states that can enter any possible decomposition of p must be of the restricted form: 

a|GHZ) +/3|W) +7 W^ , (2.81) 

with |ap + + I7P = 1. Thus, there is no need to characterize S'inv, but only to characterize 
the pure states that, under P4, are projected to p{x,y). These states are readily seen to be of the 
form: 

^/^e''^' I GHZ) + ^e''^^ |W) + 1 - x - y e'^'-" W^ . (2.82) 
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Of these, the least entangled state, for given has all coefficients non-negative (up to a global 

phase), i.e., 



\i^{x, y)) = |GHZ) + |W) + V^-x-y W) . (2.83) 
The entanglement eigenvalue of \ip{x,y)) can then be readily calculated, and one obtains 

A{x,y)= + V^t + V3(l-x-y)t2| , (2.84) 

where t is the (unique) non-negative real root of the following third-order polynomial equation: 

3y|(-t + t^) + y^{-2t^ + 1) + ^3{l- x-y){-t^ + 2t) = 0. (2.85) 

Hence, the entanglement function for \^p{x,y)), i.e., E^{x,y) = 1 — A{x,y)'^, is determined (up to 
the numerically straightforward task of root-finding). 

2.5.2 Finding the convex hull 

Recall that our aim is to determine the entanglement of the mixed state p{x,y). As we already 
know the entanglement of the corresponding pure state \^l^{x,y)), we may accomplish our aim via 
the Vollbrecht- Werner technique [2^, which gives the entanglement of p{x,y) in terms of that of 
\tp{x, y)) via the convex hull construction: Ep(x, y) = (co E^){x, y). Said in words, the entanglement 
surface z = Ep{x, y) is the convex surface constructed from the surface z = E^{x, y). 

The idea underlying the use of the convex hull is this. Due to its linearity in x and y, the state 
p(x, y) (j2.78p can [except when (x, y) lies on the boundary] be decomposed into two parts: 

p{x, y) = pp{xi,yi) + (1 - p)p(x2, 2/2), (2.86) 

with the (real, non-negative) weight p and end-points {xi,yi) and {x2,y2) related by 

pxi + {1 — p)x2 = X, (2.87a) 
pyi + {l-p)y2 = y- (2.87b) 
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Figure 2.5: Entanglement vs. the composition of the pure state \'4>{x,y)). Although not obvious 
from the plot, this entanglement surface is not convex near {x,y) = (1,0), 

Now, if it should happen that 



then the entanglement, averaged over the end-points, would give a value lower than that at the 
interior point {x,y); this conforms with the convex-hull construction. 

It should be pointed out that the convex hull should be taken with respect to parameters on 
which the density matrix depends linearly, such as x and y in the example of p{x, y). Furthermore, 
in order to obtain the convex hull of a function, one needs to know the global structure of the 
function — in the present case E^{x,y). We note that efficient numerical algorithms have been 
developed for constructing convex hulls |45j . 

As we have discussed, our route to establishing the entanglement of p{x, y) involves the analysis 
of the entanglement of \ip{x,y)), which we show in Fig. 12. 5i Although it is not obvious from the 
figure, the corresponding surface fails to be convex near the point (x,?/) = (1,0), and therefore in 
this region we must suitably convexify in order to obtain the entanglement of p{x, y). As previously 



pE^{xi,yi) + (1 - p)E^{x2,y2) < E^{x,y) 
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Figure 2.6: Entanglement of the pure state y = (1 — x)/2)) = ^fx |GHZ) + -y/(l — x)l2 |W) + 
a/(1 — x)/2 vs. X. This shows the entanglement along the diagonal boundary x + 2y = 1. Note 
the absence of convexity near x = \] this region is repeated in the inset. 



shown in Fig. 12.21 where the entanglement of \'4){x^y)) is plotted along the line {x^y) = (0, s), the 
behavior of the entanglement curve is convex. By contrast, along the line x+2y = 1 there is a region 
in which the entanglement is not convex, as Fig. 12.61 shows. The nonconvexity of the entanglement 
of \'4'{x, y)) complicates the calculation of the entanglement of p(x, y), as it necessitates a procedure 
for constructing the convex hull in the (as it happens, small) nonconvex region. Elsewhere in the 
xy plane the entanglement of p{x,y) is given directly by the entanglement of \^l'{x,y)). 

At worst, convexification would have to be undertaken numerically. However, in the present 
setting it turns out that one can determine the convex surface essentially analytically, by performing 
the necessary convexifying surgery on surface z = E^(x,y). To do this, we make use of the fact 
that if we parametrize y via (1 — x)r, i.e., we consider 



p{x, (1 - x)r) = X |GHZ)(GHZ| + (1 - x)r |W)(W| + (1 - x)(l - r)|W)(W|, 



(2.89) 



where < r < 1 [and similarly for |i/^(a;,y))] then, as a function of (x,r), the entanglement would 
be symmetric with respect to r = 1/2, as Fig. 12.71 makes evident. With this parametrization, 
the nonconvex region of the entanglement of \ip) can more clearly be identified. To convexify 
this surface we adopt the following convenient strategy. First, we reparametrize the coordinates. 
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Figure 2.7: Entanglement of the pure state (1 — x)r) ) = |GHZ) + ^/ {1 — x)r |W) + 



Y^(l — x)(l — r) vs. X and r. Note the symmetry of the surface with respect with r = 1/2. 

exchanging y by (1 — x)r. Now, owing to the Hnearity, in r at fixed x and vice versa, of the 
coefficients x, (1 — x)r and (1 — — r) in Eq. (|2.89|) . it is certainly necessary for the entanglement 
of p to be a convex function of r at fixed x and vice versa. Convexity is, however, not necessary in 
other directions in the (x, r) plane, owing to the nonlinearity of the coefficients under simultaneous 
variations of x and r. Put more simply: convexity is not necessary throughout the (x, r) plane 
because straight lines in the (x,r) plane do not correspond to straight lines in the (x,y) plane 
(except along lines parallel either to the r or the x axis). Thus, our strategy will be to convexify in 
a restricted sense: first along lines parallel to the r axis and then along lines parallel to the x axis. 
Having done this, we shall check to see that no further convexification is necessary. 

For each x, we convexify the curve z = E^{x, (1 — x)r) as a function of r, and then generate 
a new surface by allowing x to vary. More specifically, the nonconvexity in this direction has the 
form of a symmetric pair of minima located on either side of a cusp, as shown in Fig. 12.81 Thus, to 
correct for it, we simply locate the minima and connect them by a straight horizontal line. 

What remains is to consider the issue of convexity along the x (i.e., fixed r) direction for the 
surface just constructed. In this direction, nonconvexity occurs when x is, roughly speaking, greater 



46 



0.55 
. 5 

0.45 
H 0.4 

0.35 
0.3 

0.25 







0.2 



. 4 



. 6 



. 8 



Figure 2.8: Entanglement of the pure states |^(x, (1 - x)r)> = ^ |GHZ) + y/{l-x)r |W) + 
^(l-x)(l-r) W) vs. r for various values of x (from the bottom: 0.8, 0.85, 0.9, 0.92, 0.94, 0.96, 
0.98, 1). This reveals the nonconvexity in r for intermediate values of x. 
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Figure 2.9: Entanglement of the pure states |V'(x, (1 — x)r)) = ^/x |GHZ) + y^(1 — x)r |W) + 



a/(1 — x)(l — r) Wy vs. X for various values of r (from the top: 0, 0.1, 0.2, 0.3, 0.5). This reveals 
the nonconvexity in x in the (approximate) interval [0.85, 1]. 

than 0.8, as Fig. I2.9l suggests. In contrast with the case of nonconvexity at fixed r, this nonconvexity 
is due to an inflection point, at which the second derivative vanishes. To correct for it, we locate 
the point x = xq such that the tangent at x = xq is equal to that of the line between the point on 
the curve at xq and the end-point at x = 1, and connect them with a straight line. This furnishes 
us with a surface convexified with respect to x (at fixed r) and vice versa. 

Armed with this surface, we return to the (x, y) parametrization, and ask whether it is fully 
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Figure 2.10: Entanglement of the mixed state p{x,y). 



convex (i.e., convex along straight lines connecting any pair of points). Said equivalently, we ask 
whether or not any further convexification is required. Although we have not proven it, on the basis 
of extensive numerical exploration we are confident that the resulting surface is, indeed, convex. 
The resulting convex entanglement surface for p{x,y) is shown in Fig. 12.101 Figure [2.111 exemplifies 
this convexity along the line 2y + x = 1. We have observed that, for the case at hand, it is adequate 
to correct for nonconvexity only in the x direction at fixed r. 

2.5.3 Comparison with the negativity 

As introduced in Sec. II. li the negativity measure of entanglement is defined to be twice the absolute 
value of the sum of the negative eigenvalues of the partial transpose of the density matrix |14ill5l B6]. 
In the present setting, viz., the family p{x,y) of three-qubit states, the partial transpose may 
equivalently be taken with respect to any one of the three parties, owing to the invariance of p{x, y) 
under all permutations of the parties. Transposing with respect to the third party, one has 




(2.90) 



Ai<0 
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Figure 2.11: Entanglement of the mixed state y = (l-x)/2) = x |GHZ)(GHZ| + if^(|W)(W| + 
|W)(W|) vs. X. Inset: enlargement of the region x E [0.2,0.3]. This contains the only point, 
= (1/4,3/8), at which Ep{x,y) vanishes. 

where the A's are the eigenvalues of the matrix p^^, 

It is straightforward to calculate the negativity of p{x,y); the results are shown in Fig. 12.121 
Interestingly, for all allowed values of the state p{x,y) has nonzero negativity, except at 

(x,y) = (1/4,3/8), at which the calculation of the GME shows that the density matrix is indeed 
separable. One also sees from that fact that it can be obtained by applying the projection P4 to 
the (un-normalized) separable pure state ( |0) + |1) )'^^ that /9(l/4, 3/8) is a separable state. The 
fact that the only positve-partial-transpose (PPT) state is separable is the statement that there are 
no entangled PPT states (i.e., no PPT bound entangled states) within this family of three-qubit 
mixed states. The negativity surface. Fig. 12.121 is qualitatively — but not quantitatively — the same 
as that of GME. By inspecting the negativity and GME surfaces one can see that they present 
ordering difficulties. This means that one can find pairs of states p{xi,yi) and p{x2,y2) that have 
respective negativities A/i and A/'2 and GMEs Ei and E2 such that, say, A/i < A/'2 but Ei > E2. Said 
equivalently, the negativity and the GME do not necessarily agree on which of a pair of states is the 
more entangled. For two-qubit settings, such ordering difficulties do not show up for pure states, 
but can for mixed states [461 147j . On the other hand, for three qubits, such ordering difficulties 
already show up for pure states, as the following example shows: AA(GHZ) = 1 > AA(W) = 2\/2/3 
whereas for the GME the order is reversed. We note, however, that for the relative entropy of 
entanglement Eji, one has Eji{G}iZ) = log2 2 < EjiiW) = log2(9/4) [48], which for this particular 
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Figure 2.12: Negativity of the mixed state p{x,y). 
case is in accord with the GME. 

2.6 Concluding remarks 

We have considered a rather general, geometrically motivated, measure of entanglement, applicable 
to pure and mixed quantum states involving arbitrary numbers and structures of parties. In bi- 
partite settings, this approach provides an alternative — and generally inequivalent — measure to the 
entanglement of formation. For multi-partite settings, there is, to date, no explicit generalization 
of entanglement of formation. However, if such a generalization should emerge, and if it should 
be based on the convex hull construction (as it is in the bi-partite case), then one may be able to 
calculate the entanglement of formation for the families of multi-partite mixed states considered in 
the present chapter. 

As for explicit implementations, the geometric measure of entanglement yields analytic results 
in several bi-partite cases for which the entanglement of formation is already known. These cases 
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include: (i) arbitrary two-qubit mixed, (ii) generalized Werner, and (iii) isotropic states. Further- 
more, we have obtained the geometric measure of entanglement for certain multi-partite mixed 
states, such as mixtures of symmetric states. In addition, by making use of the geometric measure, 
we have addressed the entanglement of a rather general family of three-qubit mixed states analyt- 
ically (up to root-finding). This family consists of arbitrary mixtures of GHZ, W, and inverted-W 
states. To the best of our knowledge, corresponding results have not, to date, been obtained for 
other measures of entanglement, such as entanglement of formation and relative entropy of entan- 
glement. We have also obtained corresponding results for the negativity measure of entanglement. 
Among other things, we have found that there are no PPT bound entangled states within this 
general family. 

A significant issue that we have not discussed is how to use the geometric measure to provide 
a classification of entanglement of various multi-partite entangled states, even in the pure-state 
setting. For example, given a tri-partite state, is all the entanglement associated with pairs of 
parts, or is some attributable only to the system as a whole? More generally, one can envisage all 
possible partitionings of the parties, and for each, compute the geometric measure of entanglement. 
This would provide a hierarchical characterization of the entanglement of states, more refined than 
the global characterization discussed here. Another extension would involve augmenting the set of 
separable pure states with certain classes of entangled pure states, such as bi-separable entangled, 
W-type and GHZ-type states [49j. 

Although there is no generally valid analytic procedure for computing the entanglement eigen- 
value Amaxj one can give — and indeed we have given — analytical results for several elementary cases. 
Harder examples require computation, but often this is (by today's computational standards) triv- 
ial. We note that in order to find Amax for the state it is not necessary to solve the nonlinear 
eigenproblem (|2.8p : one can instead appropriately parametrize the family of separable states \(j)) 
and then directly maximize their overlap with the entangled state {ip), i.e., Amax = max^ IKc/il-f/))!!. 
Recently, Eisert and co-workers [5^ have reformulated the problem of finding the entanglement 
eigenvalue as an efficient convex optimization. Furthermore, there already exist numerical tech- 
niques for determing Ep (see, e.g., Ref. [51J), the construction of this measure being also via the 
convex-hull construction. We believe that numerical techniques for solving the geometric measure 
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of entanglement for general multi-partite mixed states can readily be developed. 

The motivation for constructing the measure discussed in the present chapter is that we wish 
to address the degree of entanglement from a geometric viewpoint, regardless of the number of 
parties. Although the construction is purely geometric, we have related this measure to entangle- 
ment witnesses, which can in principle be measured locally [H]. Moreover, the geometric measure 
of entanglement is related to the probability of preparing a single copy of a two-qubit mixed state 
from a certain pure state [33]. Yet it is still desirable to see whether, in general, this measure can be 
associated with any physical process in quantum information, as are the entanglement of formation 
and distillation. 

There are further issues that remain to be explored, such as additivity and ordering. The 
present form of entanglement for pure states, E^^^2 = 1 — A^, is not additive. However, one can 
consider a related form, i^iogj = — log2 A^, which, e.g., is additive for |V')ab l^)cD' 

^log2 ( li'l) AB ® IV'2)cD ) = ^log2 ( \i^l)AB ) + ^log^ ( \i^2)cD ) • (2-91) 

This suggests that it is more appropriate to use this logarithmic form of entanglement to discuss 
additivity issues. However, -Eiogj is not an entanglement monotone when extended to mixed states 
by convex hull, as we shall show later. 

As regards the ordering issue, we first mention a result of bi-partite entanglement measures, 
due to Virmani and Plenio [37], which states that any two measures with continuity that give the 
same value as the entanglement of formation for pure states are either "identical or induce different 
orderings in general" This result points out that different entanglement measures will inevitably 
induce different orderings if they are inequivalent. This result might still hold for multi-partite 
settings, despite their discussion being based on the existence of entanglement of formation and 
distillation, which have not been generalized to multi-partite settings. Although the geometric 
measure gives the same ordering as the entanglement of formation for two-qubit mixed states [see 
Eq. (j2.58p ]. the geometric measure will, in general, give a different ordering. It is interesting to 
note that for bi-partite systems, even though the relative entropy of entanglement coincides with 
entanglement of formation for pure states, they can give different orderings for mixed states, as 
pointed out by Verstraete and co-workers [37j . 
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We conclude by remarking that the measure discussed in the present chapter is not the same 
as the Bures measure [27j (see also Sec. II. 4p . The Bures measure of entanglement is based on the 
minimal distance between the entangled state and the set of separable mixed states. By contrast, 
the measure considered here is based upon the minimal distance between the entangled pure state 
and the set of separable pure states, and is extended to mixed states by a convex hull construction. 



53 



Chapter 3 

Connections between relative entropy 
of entanglement and geometric 
measure of entanglement 



3.1 Introduction 

In Sec. 12.11 we saw some difficulties with generalization to multi-partite settings of entanglement 
measures such as entanglement of formation and of distillation. These difficulties necessitate the 
study of other measures, e.g., the relative entropy of entanglement (REE) [27\ I48j . One of the 
reasons that the REE is interesting in the bi-partite setting is that it provides a lower bound on 
the distillable entanglement, the latter being very difficult to calculate generally. However, it is still 
non-trivial to calculate the REE for generic states, and can be more challenging for multi-partite 
states. One reason is the absence, in general, of Schmidt decompositions for multi-partite pure 
states [52j. This implies that for multi-partite pure states the entropies of the reduced density 
matrices can differ, in contrast to bi-partite pure states, as the following example shows. 

Consider a three-qubit pure state \iP)abc = a |001)+/3 |010)+7 |100), where |a|2+|/?p + |7p = 1. 
The reduced density matrices for parties A, B, and C are, respectively. 



PA 



(|a|2 + |/3p)|0)(0|+|7p|l)(l 



(3.1a) 



PB 



(I«l' + l7ni0)(0| + |/?p|l)(l 



(3.1b) 



PC 



(|7|' + |/?|2)|0)(0| + |a|2|l)(l 



(3.1c) 
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which, in general, have different entropies. Thus, for a multi-partite pure state the entropy of 
the reduced density matrix does not give a consistent entanglement measure. However, even in 
the case in which all parties have the identical entropy, e.g., a = /? = 7 = l/\/3 [53j, it is in 
general nontrivial to obtain the relative entropy of entanglement for the state. More generally, for 
pure multi-partite states, it is not yet known how to obtain their relative entropy of entanglement 
analytically. The situation is even worse for mixed multi-partite states. 

In Chapter [2] we have developed a multi-partite measure based on a simple geometric picture. 
For pure states, this geometric measure of entanglement depends on the maximal overlap between 
the entangled state and unentangled states, and is easy to compute numerically. We have examined 
this measure for several bi-partite and multi-partite pure and mixed states |37l l38] . We shall see 
in the next chapter the application to two distinct multi-partite bound entangled states [54J. In 
the present chapter, we explore connections between this measure and the relative entropy of 
entanglement. For certain pure states, some bi-partite and some multi-partite, we find that this 
lower bound is saturated, and thus their relative entropy of entanglement can be found analytically, 
in terms of their known geometric measure of entanglement. For certain mixed states, upper bounds 
on the relative entropy of entanglement are also established. Numerical evidence strongly suggests 
that these upper bounds are tight, i.e., they are actually the relative entropy of entanglement. 
These results, although not general enough to solve the problem of calculating the relative entropy 
of entanglement for arbitrary multi-partite states, may offer some insight into, and serve as a 
testbed for, future analytic progress related to the relative entropy of entanglement. 

The structure of the present chapter is as follows. In Sec. 13.21 we review the two entanglement 
measures considered in the paper: the relative entropy of entanglement and the geometric measure 
of entanglement. In Sec. 13.31 we explore connections between the two, in both pure- and mixed- 
state settings. Examples are provided in which bounds and exact values of the relative entropy of 
entanglement are obtained. In Sec. 13.41 we give some concluding remarks. The discussion in the 
present chapter is based on Ref. [55] . 



55 



3.2 Entanglement measures 

In this section we briefly review the two measures considered in the present paper: the relative 
entropy of entanglement and the geometric measure of entanglement. 

3.2.1 Relative entropy of entanglement 

The relative entropy S{p\\(t) between two states p and a is defined via 

5(p||ct) = Tr(plog2p-/>log2CT), (3.2) 

which is evidently not symmetric under exchange of p and a, and is non-negative, i.e., 5(p||(t) > 0. 
The relative REE for a mixed state p is defined to be the minimal relative entropy of p over the 
set of separable mixed states [27l [26] : 

Er{p) = min S{p\\a) = minTr (ploga p - plogg o") , (3.3) 

where P denotes the set of all separable states. 

In general, the task of finding the REE for arbitrary states p involves a minimization over all 
separable states, and this renders its computation very difficult. For bi-partite pure states, the 
REE is equal to entanglements of formation and of distillation. But, despite recent progress [56], 
for mixed states — even in the simplest setting of two qubits — no analog of Wootters' formula [22] 
for the entanglement of formation has been found. Things are even worse in multi-partite settings. 
Even for pure states, there has not been a systematic method for computing REE's. It is thus 
worthwhile to seek cases in which one can explicitly obtain an expression for the REE. A trivial 
case arises when there exists a Schmidt decomposition for a multi-partite pure state: in this case, 
the REE is the usual expression 

-J^a^logsa^, (3.4) 

i 

where the a^'s are Schmidt coefficients (with af = 1). We shall see that there exist other cases 
in which the REE can be determined analytically, even though no Schmidt decomposition exists. 
We remark that an alternative definition of the REE is to replace the set of separable states by 
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the set of postive partial transpose (PPT) states. The REE thus defined, as well as its regularized 
version, gives a tighter bound on distillable entanglement. There has been important progress 
in calculating the REE [and its regularized version, see Eq. ()3.6ip ] with respect to PPT states 
for certain bi-partite mixed states; see Refs. [57J for more detailed discussions. For multi-partite 
settings one could define the set of states to optimize over to be the set of states that are PPT with 
respect to all bi-partite partitionings. However, throughout the discussion of the present chapter, 
we shall use the first definition, i.e., optimization over the set of completely separable states. 

3.2.2 Geometric measure of entanglement 

We have introduced the geometric measure of entanglement (GME) in Sec. 12.21 in the pure-state 
settings, and generalized it to mixed-state settings via the convex- hull construction in Sec. 12. 3[ The 
essential point is to find the maximal overlap (a.k.a. entanglement eigenvalue) of the entangled 
state with unentangled states: 



where is an arbitrary unentangled pure state. The explicit form of the measure we shall be 
concerned with in the present chapter is E\og^[ilj) = — 2 log2 AmaxdV'))- We shall later show that it 
is a lower bound of the REE. 

Some of the examples we considered in previous chapter are relevant to the discussion in the 
present chapter. We now briefly recap them here. The first class contains the permutation-invariant 
pure states 



Amaxd'^)) = max|((/)|V')| 



(3.5) 




(3.6) 



which has the entanglement eigenvalue 




(3.7) 
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More generally, for n parties each a (d + l)-level system, the state 



\Sin;{k})) 



fco!fci!---fcrf! 



Permutations fcp ki 



P|0...0 ) 




(3.8) 



has the entanglement eigenvalue 



Amax("'; {k}) 




(3.9) 



The next is the totally antisymmetric state jDet^), defined via 




.in \n 



(3.10) 



which has A^^^^ = 1/n!. The generalization of the antisymmetric state to the n = p#-partite 
determinant state is 



with the 0's defined above Eq. (|238|) . The state |Det„,rf) has A^^^ = !/(#)!. 

Although the above states were discussed in terms of the GME [38], we shall, in the following 
section, show the rather surprising fact that the relative entropy of entanglement of these example 
states, is given by the corresponding expression: — 2 log2 Amax- 

3.3 Connection between the two measures 

In bi-partite systems, due to the existence of Schmidt decompositions, the relative entropy of 
entanglement of a pure state is simply the von Neumann entropy of its reduced density matrix. 
However, for multi-partite systems there is, in general, no such decomposition, and how to calculate 
the relative entropy of entanglement for an arbitrary pure state remains an open question. We now 
connect the relative entropy of entanglement to the geometric measure of entanglement for arbitrary 
pure states by giving a lower bound on the former in terms of the latter or, more specifically, via 
the entanglement eigenvalue. 
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3.3.1 Pure states: lower bound on relative entropy of entanglement 

Let us begin with the following theorem: 

Theorem 1. For any pure state with entanglement eigenvalue Ainax(^) the quantity —2 log2 Amax(V') 
is a lower bound on the relative entropy of entanglement of i.e., 

^r(|V')(V'I) > -21og2Amax(V'). (3-12) 

Proof: From the definition (j3.3p of the relative entropy of entanglement we have, for a pure state 

IV'), 

-^r(|V')(V'I) = min- (V'l logaO-IV') = -max(V'| logscrl^A) . (3.13) 

crdTJ a^T> 

Using the concavity of the log function, we have 

(V'llogaalV') <log2((Vk|V')) (3.14) 

and, furthermore, 

max log2 CT \'4}) < maxlog2((V'| IV'))) (3.15) 

although the c's maximizing the left- and right-hand sides are not necessarily identical. We then 
conclude that 

EnmiiPl) > -maxlog2((^|a|V)). (3.16) 
As any a £ V can be expanded as a = Pi | (<^i|, where |(/>i)'s are separable pure states, one has 

{^P\ a IV) = 5^P^|(</'^|V')|' < ALx(^), (3.17) 

i 

and hence we arrive at the sought result 

^^RdV'XV'l) > -21og2A^ax(V). (3.18) 



We wish to point out that such an inequality was previously established and exploited in Refs. [58J. 
When does the inequality becomes an equality? The demand that Eq. (j3.14p hold as an equality 
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implies that a (un-normalized) can be decomposed into either (a) 



a = 



(3.19a) 



where are mutually orthogonal but not orthogonal to |^), or (b) 



a=|V')(V'|+r^, 



(3.19b) 



where r (either pure or mixed) is orthogonal to ij)^ i.e., (V'l t = 0. However, the separable a 
that has either property is not necessarily the one that maximizes both sides of the inequality ()3.15p . 
unless \tjj) (and hence a) has high symmetry. On the other hand, a corollary arises from Thereom 
1 which says that for any multi-partite pure state |^), if one can find a separable mixed state a 
such that ^(pIIo-)!^^!^)^^! = -21og2 Amax( IV') ) then = -21og2 Amax( IV') )• This result follows 
directly from the fact that when the lower bound on En given in Eq. (j3.12p equals an upper bound, 
the relative entropy of entanglement is immediate. In all the examples we shall consider for which 
this lower bound is saturated, it turns out that 



is a closest separable mixed state, in which {|(/>i)} are separable pure states closest to IV")- (The 
distribution pi is uniform, and can be either discrete or continuous, and {| </>«)} are not necessarily 
mutually orthogonal.) 

We now examine several illustrative states in the light of the above corollary, thus obtain- 
ing for each of them. We begin with the permutation-invariant states \S{n,k)) of Eq. (j3.6p . 
for which A^ax was given in Eq. (13. 7p . The above theorem guarantees that EYi(y\S{n,k))^ > 
— 21og2 ^maxin, k). To find an upper bound we construct a separable mixed state 




(3.20) 




(3.21b) 



(3.21a) 
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with p chosen to maximize 1 1 (^[/^(n, /c)) 1 1 = — p)" ^, which gives p = k/n. Direct evalu- 

ation then gives 

n 

a* = ^CfcV(l-p)("-')|5(n,A;))(5(n,A:)|, (3.22) 

k=0 

and S{p\\a) = — 2 log2 Amax(f^5 fc); where p = \S{n,k)){S{n,k)\ and Ama.x{n,k) is given in Eq. (j3.7p . 
The upper and lower bounds on coincide, and hence we have that 

En{\S{n,k))) = -21og2A„,ax(n,A;). (3.23) 

The closest separable mixed state a* belongs to the case (b), i.e., Eq. (|3.19b|) . Similar equalities 
can be established for the generalized permutation-invariant n-party {d+ l)-dit states \S{n,{k})) 
of Eq. (13. Sp . We remark that the entanglements of the symmetric states \S{n,k)) (which are also 
known as Dicke states) have been analyzed via other approaches; see Ref. [39] . 

For our next example we consider the totally anti-symmetric states |Det„) of Eq. (j3.10p . It was 
shown in Ref. [40] that for these states A^^^^ = 1/n!, and hence it is straightforward to see that 
each of the n! basis states ■ ■ ■ ,in) is a closest separable pure state. Thus, one can construct a 
separable mixed state from these separable pure states [cf. Eq. (j3.20p ]: 

CJi EE ^ . . . . . . ,in\- (3.24) 

n! ^-^ 

ii,...,in 

Then, by direct calculation one gets 5'(pDot„||c"i) = ^og2{n\), which is identical to — 21og2Amax, 
as mentioned above. As in our previous examples, upper and lower bounds on coincide, and 
hence we have that £'R(|Detn)) = log2(n!). The closest separable mixed state ai belongs to the 
case (a), i.e., Eq. (j3.19ap . Similarly, for the generalized determinant state (j3.1ip one can show that 
En = log2(dP!). 

We now focus our attention on three-qubit settings. Of these, the states |S'(3,0)) = |000) 
and 15(3,3)) = |111) are not entangled and are, respectively, the components of the the 3-GHZ 
state: |GHZ) = ( |000) + |lll))/\/2. Although the GHZ state is not of the form \S{n,k)), it has 
Amax = l/\/2, and two of its closest separable pure states are |000) and |111) [38j. From these one 
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can construct a separable mixed state 



£72 = ^(|000)(000| +|111)(111|), (3.25) 



From the discussion given after Eq. (j3.18p . one concludes that £'r(GHZ) = — 21og2Amax = 1 and 
that iT2 is one of the closest separable mixed states to |GHZ). This closest separable mixed state 
(72 belongs to the case (a), i.e., Eq. (|3.19ap . With some rewriting, it can also be classified as case 
(b), i.e., 

£72 = ^|GHZ)(GHZ| + ^|GHZ-)(GHZ~|, (3.26) 

where IGHZ") = (|000) - |lll))/\/2. 
The states 



|W) = |5(3,2)) = (|001) + |010) + |100))/V3, (3.27a) 
W) = |5(3,1)) = (|110) + |101) + |011))/V3, (3.27b) 



are equally entangled, and have Amax = 2/3 [38]. Again, from the discussion after Eq. ()3.18p 
we have E'r = log2(9/4), and one of the closest separable mixed states to the W state can be 
constructed from 



^3 = I with (3.28) 



\m) ^ (V273|0) + e**^/l73|l))^^ (3.29) 

which gives the result 

a, = ^|W)(W| + ^|W)(W| + ^|000)(000| + 1|111)(111|. (3.30) 

We remark that the mixed state cr^ is not the only closest separable mixed state to the W state; 
the following state C74 is another example (as would be any mixture of and a^): 

2 

C74 ^ ^ ^ m2nk/3)){iji2nk/3)\ = ||W)(W| + ^|W)(W| + h^iCl (3.31a) 



3 

fc=0 
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Figure 3.1: The solid curve represents Eiog^{s) of the pure state ^/s\W) + y/1 — s |W) vs. s. The 
dots are the corresponding relative entropies of entanglement, obtained numericahy. 

where 3 |^) = 2\/2 |000) + |111). These closest separable mixed states of W state belong to the case 
(b), i.e., Eq. (I3l9bll . 

Having obtained REE for W and W, it is interesting to examine the REE of the following 

WW(s) ) = ^/s |W) + a/1 - s W). We have not been able to find an 



superposition of the two: WW(s)^ = ^/s |W) + a/1 — s ^ 
analytical result for REE, but we can compare the analytical expression for — 21og2 Amax(WW(s)) 
with the numerical evaluation of £'r(WW(s)), and we do this in Fig. 13. li As we see in this figure, 
the qualitative behavior of the two functions is similar, but — 21og2 A^ax and £'r only coincide at 
the two end-points, s = and s = 1. 



3.3.2 Mixed states: upper bound on relative entropy of entanglement 

In Chapter [2] the procedure was given to find the geometric measure of entanglement, E^^^2, for the 
mixed state comprising symmetric states: 



p{{p}) = ^Pk\S{n,k)){S{n,k)\. 



(3.32) 



Here, we focus instead on the quantity -E^iogj ; but the basic procedure is the same. The first step is 
to find the entanglement eigenvalue A„({g}) for the pure state 



\S{n,k)) , 



(3.33) 
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thus arriving at the quantity 

f({g}) = -21og2An(W). (3.34) 
Then the quantity -Eiogj for the mixed state (13.321) is actually the convex hull of the expression ()3.34p : 

^iog,(p(M))=cof(M). (3.35) 

This prompts us to ask the question: Can we find REE for the mixture of \S{n, k)) in Eq. (j3.32p ? 
To answer it, we shall first construct an upper bound to REE, and then compare this bound with 
the numerically evaluated REE. To accomplish the first step, bearing in mind the fact that any 
separable mixed state will yield an upper bound, we consider the state formed by mixing the 
separable pure states ICiOA)) [cf- Eq. IK22t\ ]: 

^ie)= [ ^|C(^,</>))(e(^,0)|=X;CfcCos2'=^sin2("-'^)0|S(n,fc))(5(n,A;)|, (3.36) 

where 

\^{9,4>)) = (^cos0|O) +e''^sin0|l))®". (3.37) 
We then minimize the relative entropy between p{{p}) and cr{6), 

S{p{{p})\Wm=j:P'^log J\ (3.38) 

^ C J: cos^*^ V sm ' 



with respect to 9, obtaining the stationarity condition 



tan^^^^^#^. (3.39) 



Due to the convexity of the relative entropy. 



S ^QiPiW'^qm < ^qiS{pi\\ai), (3.40) 

\ i i / i 

we can further tighten the expression of the relative entropy by taking its convex hull. (Via the 
convexification process, i.e., the convex hull construction, the corresponding separable state can 
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also be obtained.) Therefore, we arrive at an upper bound for the relative entropy of entanglement 
of the mixed state p{{p}): 

Ej,{p{{p}))<coF{{p}), (3.41) 

where 

F(M) = ^p, log, = E^^- Cj:aHn-ar-^^ ^''''^ 

where the angle 6 satisfies Eq. (|3.39p . = n\j {k\{n — /c)!), and a = J2kPk ^■ 

Having established an upper bound for REE for the state p{{p}), we now make the restriction 
to mixtures of two distinct n-qubit states |5(n,A;i)) and |S'(n, A;,)) (with ki / A;,): 

Pn-MMis) = s\S{n, h)){S{n, A:i)| + (1 - s)\S{n, k2)){S{n, A:,)!- (3.43) 

One trivial example is Pn;0,n{s), which is obviously unentangled as it is the mixture of two separable 
pure states |0®") and ll*^"). Other mixtures are generally entangled, except possibly at the end- 
points s = or s = 1 when the mixture contains either \S{n,0)) or \S{n,n)). We first investigate 
the two-qubit (i.e. n = 2) case. Besides the trivial mixture, P2;0,2, there is only one inequivalent 
mixture, P2;0,i(s) [which is equivalent to /92;2,i(s)], which is — up to local basis change — the so-called 
maximally entangled mixed state |59l [36] (for a certain range of s) 

P2;o,i = s |11)(11| + (1 - s)|^+)(^+|, (3.44) 

where |^'+) = (|01) + 10)/V^. The function F for this state [denoted by -F2;0,i(s)] is 

i^2;0,i(s) = s log2 (Y^^ + (1 - s) loga , (3.45) 

which is convex in s. It is exactly the expression for the relative entropy of entanglement for the 
state P2;0,i found by Vedral and Plenio [26] (see their Eq. (56) with A replaced by 1 — s). 

For n = 3 there are three other inequivalent mixtures: P3-Q^i{s) [equivalent to P3;3,2('S)]! P3;0,2{s) 
[to /03;3,i(s)], and P3-i^2is) [to P3;2,iis)]- In Fig. 13.21 we compare the function F in Eq. (|3.42p . 
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its convex hull coF, and numerical values of En obtained using the general scheme described in 
Ref. |i26] extended beyond the two-qubit case. The agreement between co-F and the numerical 
values of £'r, appears to be exact. 

For n = 4 there are five inequivalent nontrivial mixtures: /94;o,i(s), /04;0,2('S), /O4;0,3(s)) P4;i,2is), 
and P4;i,3(s). In Figs. 13.31 and 13.41 we again compare the function F in Eq. (j3.42p . its convex hull 
coF, and numerical values of E'r. Again the agreement between coF and the numerical values of 
Er appears to be exact. 

From these agreements, we are led to the following conjecture: 
Conjecture 1: The relative entropy of entanglement Er {p{{p})) for the mixed states p{{p}) is given 
exactly by coF{{p}). 

For the states that we have just considered, we now pause to give the formulas for Er suggested 
by the conjecture. For the three-qubit mixed state P3;2,i('S), its conjectured Er is 

9s 9(1 — s) 

^^Qg2 . ^ + (1 - s)log2 7;^ TITT^T- (3.46a) 

{l + sy{2-s) (2-s)^(l + s) 



For P3;0,i{s), it is 



For P4;0,i(s), it is 



For P4;i,2{s), it is 



For P4;i,3{s), it is 



27s 9 

^^°^^(2T^ + (^-^)^°^^(2T^- ^'-'^^^ 



256s 64 
s^og2 + 1-s log2 3 - (3.47a 

(3 + s)* (3 + s)'^ 



64s . 128(1 -s) 

(2-s)(2+s)3 + (1-^) 3(2_,)2(2+,)2- (3-47b) 



64s 64(1 — s) 

(3-2s)(l + 2s)3 + (1-^)^°^^ (3-2s)3(l + 2s) - ^^-'"'^ 



For states such as P3;o,2) Pi;0,2, and P4;o,3, convexifications (i.e. convex hull constructions) are 
needed; see Figs. 13.21 13.31 and 13. 4[ In Fig. 13.51 we give an example of a seven-qubit state, viz., 

P7;2,5(s)- 

Although we have not been able to prove our conjecture, we have observed some supporting 
evidence, in addition to the numerical evidence presented above. We begin by noting that the 
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0.2 0.4 0.6 0.8 1 



Figure 3.2: Comparison of F (solid curve), coF (convexification, indicated by dashed line) and the 
numerical value of En (dots) for the states /53;0,i(s), /93;0,2(s), and P3;i,2(5) (from top to bottom). 
Note that the log function is implicitly base-2. 
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Figure 3.3: Comparison of F (solid curve), its convex hull (dash line), and the numerical value of 
E^i for the state P4;0,3(s). Upper panel shows the whole range s G [0, 1], whereas the lower panel 
shows a blow-up of the range s G [0, 0.01]. 
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Figure 3.5: The function F (solid curve) and its convex hull (dashed line indicates convexification) 
for the seven-qubit mixed state p7-2,5is). 

states p{{p}) are invariant under the projection 



with U{(t>){ |0) , |1) } ^ { |0) ,e~'^ |1) }. Vollbrecht and Werner [24j have shown that in order to 
find the closest separable mixed state for a state that is invariant under projections such as P, it 
is only necessary to search within the separable states that are also invariant under the projection. 
We can further reduce the set of separable states to be searched by invoking another symmetry 
property possessed by p{{p})- these states are also, by construction, invariant under permutations 
of all parties. Let us denote by Ilj one of the permutations of parties, and by nj(/j) the state 
obtained from p by permuting the parties under Ilj. We now show that the set of separable states 
to be searched can be reduced to the separable states that are invariant under the permutations. 
To see this, suppose that p is a mixed state in the family ()3.32p . and that a* is one of the closest 
separable states to p, i.e., 




(3.48) 



En{p)^minS{p\\a) = S{p\\a*). 



(3.49) 



As p is invariant under all Ilj, we have 




(3.50) 
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where A^n is the number of permutations. By using the convexity of the relative entropy we have 

Er{p) > S (^p\\[Y.^^icTl/Nn]^ . (3.51) 

However, because of the extremal property, Eq. ()3.49p . the inequality must be saturated, as the 
left-hand side is already minimal. This shows that 

^"^]^En^(^*) (3.52) 

i 

also a closest separable mixed state to p, and is manifestly invariant under all permutations. Thus, 
we only need to search within this restricted family of separable states. 

It is not difficult to see that the set T>s of all separable mixed states that are diagonal in the 
basis of {\S{n,k))} can be constructed from a convex mixture of separable states in Eq. (j3.36p . 
That is, for any fx^ G Ds we have a decomposition 

a, = ^Uai9^), (3.53) 

i 

where ti > 0, J2i'^i — ^^'^ cr(^j) is of the form (j3.36p . This is because the separability of the 
states (j3.32p implies that there exists a decomposition into pure states such that each pure state 
is a separable state. Furthermore, because {\S{n,k))} are eigenstates of p{{p}), the most general 
form of the pure state in its decomposition is 

Y,VQ^e"^^\Sin,k)). (3.54) 

k 

This pure state is separable if and only if it is of the form (I3.37P , up to an overall irrelevant phase. 
As p{{p}) is invariant under the projection P (|3.48p . a pure state in Eq. (j3.37p will be projected to 
the mixed state in Eq. (j3.36p under P. Thus, every separable state that is diagonal in {\S{n,k))} 
basis can be expressed in the form (j3.53p . 

Hence, our conjecture (j3.4ip ensures (via any necessary convexification) that it is at least the 
minimum (of the relative entropy) when the separable mixed states are restricted to Ds- However, 
in order to prove the conjecture, one would still need to show that the expression is also the 
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Figure 3.6: Comparision of E (dashed curve) and F (solid curve) for the eleven-qubit mixed state 
Pii;2,6('S)- 

minimum when the restriction to is relaxed. 

We remark that our conjecture is consistent with the results of Ishizaka [60], in that our conjec- 
tured a* satisfies the condition that [p, o"*] = and that g* has the same reduction as p for every 
party. Furthermore, suppose a* (diagonal in the basis {|5'(n, A;))}) represents the separable state 
that gives the conjectured value of REE: 

a* = Y,rk\S(n,k)){S(n,k)\, (3.55) 

k 

where the r's can be obtained by finding the convex hull of the function F in Eq. (j3.42p . Now 
consider any separable state r in the Hilbert space orthogonal to the subspace spanned by {S{n, k)}. 
We need to show that the separable state a{x) = xa* + (1 — x)t, for any x £ [0, 1], gives greater 
relative entropy with p{{p}) in Eq. (|3.32p than a* does with p{{p}), i.e., 

S ip{{p})\\a{x)) > S {p{{p})\\a*) . (3.56) 




Writing out the expression explicitly, we have that 

S (p(W)lla(x)) = J^p, log ^ > ^pfc log ^ = 5 (p(M)lk*) . (3.57) 

^ — ^ X ri. ^ — ' ri. 
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Note that r gives no contribution in the relative entropy, as it is orthogonal to p{{p}), and that 
we have not used the fact that r is separable. But to prove Conjecture 1 we need to show that 
Eq. (j3.56p holds if separable r is not orthogonal to the subspace spanned by {S{n, k)}. 

Recall that for pure states we found the inequality -Eiogj ^ -E'r- Does this inequality hold for 
mixed states? We do not know the complete answer to this question, but for the mixed state p{{p}) 
we shall at least find that this inequality would hold if Conjecture 1 holds. To see this, we first 
establish that £{{q}) is a lower bound on F({q}); see the example in Fig. 13. 6[ The proof is as 
follows. Recall that 



= -21og2 



max y/pj; cos^' 6 sin" ''^ i 



(3.58a) 



By the concavity of log, we then have 



-21og2 



E 



Pk JC^ cos^ e sin'^-^ i 



< 



Pk 



C^cos2^'esin2("-^) 



(3.58b) 



Hence 



min —2 log2 



E 



SyCf cos^'^sin"-^ 



<min^p,,log2 — 



Pk 



C^cos2^esin2("-'=) 



-, (3.58c) 



or equivalently 



£i{p}) < Fi{p}). 



(3.58d) 



If Conjecture 1 is correct then by taking the convex hull of both sides of this inequality we would 
have 



(3.58e) 



for the family of states (j3.32p . Notice that we have also shown that this relation holds for arbitrary 
pure states. It would be interesting to know whether it also holds for arbitrary mixed states. 
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3.4 Concluding remarks 

We have provided a lower bound on the relative entropy of entanglement for arbitrary multi-partite 
pure states in terms of their geometric measure of entanglement. For several families of pure states 
we have shown that the bound is in fact saturated, and thus provides the exact value of the relative 
entropy of entanglement. For mixtures of certain permutation-invariant states we have conjectured 
analytic expressions for the relative entropy of entanglement. 

It is possible that our results on the relative entropy of entanglement might be applicable to 
the checking of the consistency of some equalities and inequalities |48t [611 [62] regarding minimal 
reversible entanglement generating sets (MREGSs). Consider, e.g., the particular family of n-qubit 
pure states {\S{n,k))}, the relative entropy of entanglement of which we have given in Eq. (|3.23p . 
Now, if we trace over one party we get a mixed {n — l)-qubit state: 

TTi\S(n,k)){S(n,k)\ = '^\S(n-l,k)){S(n-l,k)\ + -\S(n-l,k-l)){S(n-l,k-l)\. (3.59) 

n n 

We have also given a conjecture for the relative entropy of entanglement for this mixed state. If 
we trace over m parties, the reduced mixed state would be a mixture of — m, g))} [with 

q < {n — m)], and again we have given a conjecture for its relative entropy of entanglement. For 
example, if we start with |5'(4, 1)), and trace over one party and then another, we get the sequence: 

|5(4, 1)) ^ /93;0,l(l/4) - P2;0,l(l/2), (3.60) 

for which we have given the corresponding relative entropies of entanglement in Eqs. (|3.23|) . (|3.46bp 
and (|3.45p . (To be precise, the second formula is a conjecture; the others are proven.) The afore- 
mentioned equalities and inequalities concerning MREGS usually involve only the von Neumann 
entropy and the regularized (i.e. asymptotic) relative entropy of entanglement of the pure state and 
its reduced density matrices. The regularized relative entropy of entanglement is defined as 

E^ip)^ hm -i?R(p«-). (3.61) 

n^oo n 

The calculation of the regularized relative entropy of entanglement is, in general, much more difficult 
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than for the non-regularized case, and the (in)equalities involving the regularized relative entropy 
of entanglement are thus difficult to check. Nevertheless, it is known that E^^ < Er, so we can 
check their weaker forms by replacing E^ by Er, and the corresponding (in) equalities by weaker 
inequalities. 

Plenio and Vedral |48j have derived a lower bound on the REE of a tripartite pure state pabc = 
|^)(^| in terms of the the entropies and REE's of the reduced states of two parties: 

max{£;R(pAB) + S{pab), Er{p Ac) + S{pAc), Er{pbc) + S{pBc)} < Er{pabc), (3.62) 

where pab = Trc(/9ABc) (and similarly for pAC and PBc) and S{p) = — Trplog2/9 is the von 
Neumann entropy. They have further found that this lower bound is saturated by |GHZ) and |W). 
This raises an interesting question: is the above lower bound (for n-partite pure states) saturated 
by the states that saturate the lower bound -Eiogj = — 2 log2 Amax(V') ^ Er{iP)7 Numerical tests 
seem to suggest that the Plenio- Vedral bound is tighter than -Eiogj- If this is the case then all 
states that saturate the lower bound -Eiogj oii will saturate the Plenio- Vedral bound. Based on 
Conjecture 1, we can show that for p\2...n — 

\S{n,k)){S{n,k)\ the inequality 

max{£;R(p^2 i n) + 'S'(Pi2 I „)} < £^R(Pi2...n) (3.63) 

is saturated, where p^^ 1 n = Trj(pi2...n) is the reduced density matrix obtained from pi2...n by 
tracing out the i-th party. The proof is as follows. As \S{n,k)) is permutation- invariant, there is 
no need to maximize over all parties, and we can simply take i = 1, obtaining the reduced state 
Pn-i;k-i,k{k / n) as in Eq. (j3.59p . As the corresponding function Fn-i-k-i,k{s) of Pn~i;k-i,k{s) is 
convex for s G [0,1], we immediately obtain from Conjecture 1 that, for Pn~i;k~i,k{k/n), 



ER{pjj^i.k-i,k{k/n)) = log2 



k\ / n—k 



k , k n — k ^ n — k ^, . 

+ - log2 - H log2 (3.64a) 

n n n n 



n J \ n 

= ER{\S{n,k)))-S{p^i,k-i,k{k/n)). (3.64b) 

Therefore, the bound in Eq. ()3.63p is saturated for pi2...n = \S{n,k)) {S{n,k)\. 

A major challenge is to extend the ideas contained in the present Paper from the relative 
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entropy of entanglement to its regularized version, the latter in fact giving tighter upper bound on 
the entanglement of distillation than the former in the bi-partite settings. The alternative way of 
defining the relative entropy via the optimization over PPT states may also been used, in view of 
the recent progress on the bi-partite regularized relative entropy of entanglement [57] . 

We now explore the possibility that the geometric measures can provide lower bounds on yet 
another entanglement measure — the entanglement of formation. If the relationship E'r < Ep 
between the two measures of entanglement — the relative entropy of entanglement E'r and the 
entanglement of formation Ep — should continue to hold for multi-partite states (at least for pure 
states), and if Ep should remain a convex hull construction for mixed states, then we would be able 
to construct a lower bound on the entanglement of formation: 

Eiog^ip) = min Vpi^^iogjdV'j)) < min Vpi^R(IV'i)) 

< mmy^p,Ep{\i;i)) = Ep{p), (3.65) 

where {pi} and {ijji} are such that p = Yl,iP'i\'^i){'^i\- Thus, E\o^^{p) is a lower bound on Ep{p). 
By using the inequality (1 — x^)log2e < — 21og2X (for < x < 1), one further has has that 

{\og^e)E^Ap)<ExogM<Mp)- 

We remark that E^^^2 has been shown to be an entanglement monotone [371 138] . i.e., it is not 
increasing under local operations and classical communication (LOCC). However, E\o^^ is not a 
monotone, as the following example shows. Consider the bi-partite pure state 

l« ^ 7m l°"> + Tfhw^ 111) + 122) + ... + m) ), (3.66) 

with l^l < 1, for which -Eiogj = log2(l + Nx^). Suppose that one party makes the following 
measurement: 

7Wi = |0)(0|, 7^2 = |l)(l| + |2)(2| + --- + |iV)(iV|. (3.67) 

With probability Pi = 1/(1 + Nx'^) the output state becomes IV'i) = |00); with probability P2 = 
Nx'^/{1 + Nx'^) the output state becomes IV'2) = + |22) + ••• + |iViV))///V, for which 
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Figure 3.7: The function f{4,x). It shows the violation of monotone condition (j3.69p when the 
function is negative. 

-E'logj = log2 N. For -Eiogj to be a monotone it would be necessary that 

Eiog.W > Pi^iog^lV'i) + P2Eiog,{i^2)- (3.68) 

Putting in the corresponding values for the P's and E'logj's, we find that this inequality is equivalent 
to 

/(iV, x) ^ log2(l + iVx^) - log2 iV > 0. (3.69) 

As this is violated for certain values of x with > 2, as exemplified in Fig. 13.71 for the plot of 
/(4, x), we arrive at the conclusion that -Eiogj is, in general, not a monotone. 

We conclude by mentioning that certain results reported in the present chapter have recently 
been applied by Vedral [B^ to the macroscopic entanglement of ry-paired superconductivity. 
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Chapter 4 

Bound entanglement 



4.1 Introduction 

As we have discussed above, we arc motivated to study the quantification of entanglement for 
the basic reason that entanglement has been identified as a resource central to much of quantum 
information processing. As we have also discussed above, to date, progress in the quantification of 
entanglement for mixed states has resided primarily in the domain of bi-partite systems. For multi- 
partite systems in pure and mixed states, the characterization and quantification of entanglement 
presents even greater challenges. 

We have introduced previously the notion of the entanglement of distillation, entanglement cost, 
and entanglement of formation. The entanglement of distillation {E-q) is the optimal asymptotic 
yield of Bell states, given an infinite supply of replicas of an identical quantum state shared between 
two distant parties. The entanglement cost {Eq) is the minimum asymptotic cost ratio of Bell 
states consumed to create an large number of copies of a certain quantum state shared between two 
distant parties. For pure states shared between two distant parties the two processes are reversible, 
i.e., £'d = Ec- For mixed states it is expected that Ed < Eq; otherwise recycling the two 
processes would churn out more entanglement than there was initially. The quantity entanglement 
of formation Ep is an averaged version of Eq- 

All these three quantities, £'d, Eq, and Ep, especially the first two, are very difficult to calculate. 
The calculation of these quantities for general quantum states remains a challenge in quantum 
information theory. Another challenge regarding Ep concerns the question of whether or not Ep, 
defined in an average sense, is equal to Ec, which is defined for asymptotically large numbers of 
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copies. This is the so-called additivity problem for entanglement of formation, mentioned earlier. 

There is, however, some progress in small dimensions, especially for two-qubit systems. Woot- 
ters' formula |22j for two-qubit entanglement of formation is the most prominent example. It sets 
an upper bound on the entanglement cost. It is obvious that an entangled state has nonzero Eq 
and Ep. It was also shown by Horodecki and co-workers |18j that all two-qubit entangled states 
can be distilled into Bell states, and hence have nonzero ii^D- It was then thought that all entan- 
gled states, however small the entanglement, could be distilled. But shortly after, Horodecki and 
co-workers \19\ 164] showed that in dimensions higher than two-qubit (C^ C"^) and qubit-qubtrit 
{C"^ ® C^) [such as and ® C^] there exist entangled states that have zero entanglement 
of distillation. These states, however entangled, cannot be distilled into any pure entangled states. 
The entanglement used to create them is somehow bound and inextractable! Such states are called 
hound entangled states, and their entanglement is called bound entanglement. The bound entan- 
gled states that Horodecki constructed are bi-partite and have zero negativity, i.e., have positve 
partial transpose (PPT). Bound entanglement is not limited to bi-partite states. Bennett and co- 
workers |65) constructed both bi- and multi-partite bound entangled states using the technique of 
unextendible product bases; see also Appendix iFl Other multi-partite bound entangled states have 
also been found [66j, including the two ^GTIEHIIMI that we shall discuss in this chapter. 

Although bound entanglement seems useless, as it cannot be used alone for quantum commu- 
nication, Horodecki and co-workers [70] found that, surprisingly, a certain threshold fidelity of a 
teleportation process that cannot be achieved via a single copy of non-maximally entangled pure 
state can actually be achieved if combined with a supply of bound entanglement. Recently, Shor 
and co-workers [68] have shown that two multi-partite bound entangled states (of Smolin's), ten- 
sored together, can be be distilled (via LOCC) into a Bell state shared between two of the parties. 
Bound entanglement clearly has richer properties than was thought initially. 

We have seen that bound entangled states are states that are entangled, but from which no 
pure entangled state can be distilled, provided all parties are allowed only local operations and 
classical communication (LOCC). The distillable entanglement (-Ed) is thus zero by definition. 
Bound entangled states can be either bi-partite or multi-partite, the latter possibly exhibiting more 
structure than the former. However, it does take nonzero entanglement to create bound entangled 
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states under LOCC. But exactly how much entanglement is bound in these strange states is not 
known analytically. In the present chapter, we study the entanglement content of two distinct 
types of bound entangled state: Smolin's four-party unlockable bound entangled state [67\ [68] 
and Diir's A^-party Bell-inequality- violating bound entangled states [69]. For each, we determine 




shown in the previous chapter that, under certain circumstances, these give lower bounds on their 
multi-partite Ep. In particular, we showed that i?gi^2 (p) log2 e < E\o^^[p) < Ep{p). In addition, 
we make conjectures concerning the relative entropies of entanglement for these bound entangled 
states. Although quantities such as the geometric measure or the relative entropy of entanglement 
may not be able to reveal the exact nature of bound entanglement, they nevertheless quantify for 
these bound entangled states the content of entanglement that is inextractable. The discussion in 
the present chapter is based on Ref. [5^ . 

We now turn to the calculations of entanglement for the two bound entangled states: Smolin's 
and Diir's. 

4.2 Smolin's four-party unlockable bound entangled state 

Consider the four-qubit mixed state 



where the |^)'s are the four Bell states: (|00) ± |ll))/\/2 and (|01) ± |10))/\/2. Now, the state 




analytically their geometric measure of entanglement E^^^i and the related quantity Ei^g^ . We have 




(4.1) 



,ABCD 



can be conveniently rewritten as 




(4.2) 



j=0 



where |X)'s are the four orthogonal GHZ-like states: 



Xo) = -^ (10000) + 11111)), |Xi) 
X2) = ^(|0101) + |1010)), 1x3) 



V2 



V2 



1 



1 



(|0011) + |1100)) 



(|0110) + |1001)). 
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From the decomposition in Eq. (j4.2p we immediately see that the state 

pABCD ig 

invariant under 

any permutations of the four parties. 

If any two of the four parties, say C and D, get together, they can do a Bell- measurement 
(namely, measurement done in the basis of |^'j)'s). Depending on the result i = 0, ...,3, they 
broadcast the outcome to A and B, from which A and B can then establish a definite Bell state. 
This shows that the state p^^'^'^ must be entangled. But if all the four parties are far apart, 
they have no way of distilling any pure entangled states. This can be seen from the fact proven 
in Ref. [20] that if two parties are on opposite sides of a separable cut, then they will remain in 
a separable cut under any local quantum operations and classical communication. From the form 
in Eq. (14. ip we see that A cannot establish any entanglement between either C or D, as A is in 
the opposite side of a separable cut from C and D. But from Eq. ()4.2p we know that the state is 
invariant under any permutation of the parties, hence, A cannot establish any entanglement with 
B (by exchanging C with B), either. Therefore, the state 

pABCD 

bound entangled. 

Our goal here is to calculate how much entanglement is bound in the state. As the state is 
bound entangled, it is equivalent to ask how entangled is the state? For the purpose of using GME 
to quantify entanglement, we need to characterize all decompositions of the mixed state into pure 
states. The most general decomposition of a mixed state p into pure states can be expressed as 

M n 
P = '^\'fk){'fk\, with \(pk) = y^Uki\f>^i , (4.3) 

k=l i=l 

where is an integer not smaller than n, the number of orthonormal eigenvectors {l^^j)} (with 
nonzero eigenvalues {Aj}) of p, |(^)'s are un-normalized, and Li satisfies X^^i^fci^fej = ^ij- Thus, 
the most general pure state that appears in the decomposition of Smolin's state is 

\'Pk) = ^^Uki\Xi). (4.4) 

1=0 

Our goal is to minimize YlkPk ^puve{\^k)) over all possible Ws, where -Bpure is some pure-state 
entanglement {E^^^2 or -Eiogj in onv considerations), pk = {(pk\(pk)-, and \(pk) is the normalized state 
Wk) = \'^k) I \fPk- Making a general minimization for an arbitrary mixed state is extremely difficult. 
However, for the mixed state p^^'^^ we shall show that the decomposition in Eq. (j4.2p does indeed 
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minimize the average entanglement over pure-state decompositions. As in Eq. ()4.4p . \ip) can be 
explicitly written as \ip) = Yli=oV^^^'^^ l^j)> where the q^s are non-negative, satisfying = 1, 

and the i;^'s are phases. For fixed g's, the state has a maximal entanglement eigenvalue when all 
phases are zero. We shall show shortly that its maximal entanglement eigenvalue is l/\/2i which 
is achieved by the |^)'s. 

The entanglement eigenvalue of the state \ip) = X]?=o maximal overlap with the 

separable state |$) = (^j^-^(^Ci |0) -|- Si |1)), where Cj = cos0j and Sj = sin6'j with < < 7r/2. 
Thus 

{^\ip) = ^/ qo/2 (C1C2C3C4 -I- S1S2S3S4) + V 91/2 (C1C2S3S4 + S1S2C3C4) 
+ ^^2/2 (C1S2C3S4 + S1C2S3C4) + \/93/2 (C1S2S3C4 + S1C2C3S4), 

which has maximum I/V2. To see this, use the Cauchy-Schwarz inequality, treating as one vector 



{\/9o/2,\/gi/2, Y^g2/2,Y^<?3/2| (whose modulus is l/\/2), and the corresponding coefficients as 
another vector, whose modulus can be shown to be no greater than 1: 

(C1C2C3C4 + 81828334)^ + (C1C2S3S4 + SiS2C3C4)^ 

+ (C1S2C3S4 + 81028384^ + {C18283C4 + 81020384)^ < 1. 

By subtracting the left-hand side from 1 and making some algebraic manipulation, we arrive at the 
non- negative expression (hence the sought result): 

(C1C2C3S4 - 81828364^ + (C1C2S3S4 - 81820384^ + 
(C1S2C3C4 - 81028384^ + (S1C2C3C4 - 01828384^ > 0. 



The states \Xi) are GHZ-like states and have Amax = l/v2 and they clearly saturate < 
I/V2. Hence, we have 

^.n^(/^^^) = l/2, E,^,^{p^^^^) = l. (4.5) 
This suggests that although bound entangled, Smolin's state has a very high degree of entanglement. 
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the same as that of a 4-partite GHZ state. This high degree of entanglement seems to manifest in 
some bi-partite partitioning, e.g., {A:BCD} (as we discuss below). 

To compare the results with other measures of entanglement, we conjecture (and later prove) 
that Eji = 1 for this state and one of its closest separable mixed states is 

^(|0000)(0000| + |1111)(1111| + |0011)(0011| + |1100)(1100| 
8 

+|0101)(0101| + |1010)(1010| + |0110)(0110| + |1001)(1001|). 

The negativity M (a value used to quantify the degree of bi-partite inseparability of states and 
defined as twice the absolute sum of negative eigenvalues of the partial transpose (PT) of the density 
matrix with respect to some bi-partite partitioning) is zero for any 2/2 partitioning, e.g., {AB : 
CD}, but nonzero for 1/3 partitioning, e.g.,{A:BCD}. (This nonzero negativity also demonstrates 
that the state entangled.) Specifically, A/aiBCD = 1 but A/'ab:CD = 0. 

Let us now turn to Diir's bound entangled states. 



4.3 Diir's A^-party bound entangled states 

Diir j69j found that for iV > 4 the following state is bound entangled: 

PN^J^(^^G){^G\ + lYl{Pk + Pk)^ , (4.6) 

where = ( [O'^^) + e'°'^ is a iV-partite GHZ state; Pk = \uk){uk\ is a pro- 

jector onto the state \uk) = |0)-^ [0)2 • . . . . . |0)^; and Pk = \vk){''^k\ projects onto 1?;^) = 
11)2 . . . |0);, . . . For > 8 this state violates the Mermin-Klyshko-Bell inequality |69] : 

violation was pushed down to A > 7 by Kaszlikowski and co-workers [7l] for a three-setting Bell 
inequality; it was pushed further down to A > 6 by Sen and co-workers |72j for a functional Bell 
inequality. For these inequalities, see Appendix [El for more detail. These results are interesting and 
somewhat surprising, as one might expect that bound entangled states has low entanglement that 
they could not violate any Bell inequality. So how entangled are Diir's bound entangled states? 
We remark that the phase a^v in \^g) can be eliminated by local unitary transformations, and 
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hence we shall take ckjv = in the following discussion. In fact, if we consider the family of states 

1 - ^ 

Pm{x) = x\^g){^g\ + ^ E (^*^ + ^k), (4.7) 

fe=i 

we find that for N > 4 the state is bound entangled if < x < 1/{N + 1) and is still entangled 
but not bound entangled if x > 1/{N + 1). This can be seen from the fact that the negativities of 
Pn{x) with respect to the two different partitions (1 : 2 • • • AT) and (12 : 3 • • • iV) are 

Afi:2-N{pN{x)) =ma^{0,[{N+l)x-l]/N}, (4.8a) 
Mi2:3-n{pn{x)) = X. (4.8b) 

By applying arguments similar to those used to calculate entanglement for Smolin's state, we 
have that the general pure state in the decomposition of Pn{x) is 

N 

^e"^° \^g) + (V^^'^' 1"^) + Vne"*''^ \v^)), 

k=l 

where q's and r's are non-negative and satisfy ^^iqk + fk) = 1- In this family, the state with the 
least entanglement (or maximum Amax) for fixed {y,qk,rk} is the one with all phase factors zero: 

N 

I* (y, U, r})) = ^/y \^g) + 0^ 1"^) + l^i) ) • 

k=l 

Next, we ask: For fixed y, what is the least entanglement that the above state can have? Take a 
separable state of the form |$) = ®^^i{ci |0) + |1) ); its overlap with |^'(y, {q,?"})) is then 

N 

(*|$) = \/yj2{ci ■ ■ ■ CN + Si - ■ ■ Sn) + \/l-yy~^(-\/gfcCl • • • Sfe • • - CAT + V^Sl • • - Cjfc • • • sn)- 

k=l 



This can be shown to no greater than (2 — y)/2, again by a Cauchy-Schwarz inequality, taking 

{V^,Wi^ - y)qk},Wi^ - y)rk}} 

as the first (2A^-|-l)-component vector (with modulus -^/(2^^y)72) and the corresponding coefficients 
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as the second one, whose modulus can be shown to be no greater than 1 for > 4: 

TV 

fN= {ci - ■ -CN + Si - ■ ■ Sn) + ^ {(Cl • • • Sfc • • • Cat)^ + (Sl • • - Cfc • • • Sn)"^} < 1- 

k=l 

First, making similar arguments as previously, one can show that /4 < 1. One can also show that 
In+i < In- Thus by induction, we have proved the inequality. 
The bound can be saturated, e.g., by 

\^±,u,k{y)) = VvI^g) ±V^^\uk) , (4.9a) 

\^P±,v,kiy)) = Vy\'^G)±V^\vk), (4.9b) 

for which Aiiiax(y) = y^(2 — y)/2 [73]. This can be seen as follows. As one can make local relative 
phase shifts to transform ^{"^g) + \/l — y \uk) to y^l^'c) — ~ V \'^k)-, they have the same 
entanglement. The change from ^\^g) =t ^1 — y \uk) to ^\^g) =t ^f^^--y\vk) is simply a flipping 
of to 1, and vice versa. The mapping from k to k' is just a relabelling of parties. Thus, we need 
only consider the state 

v^(ioo---o) + |ii---i)) + yr^iio---o). 

As this state is invariant under permutation of all parties except the first one, and as the coefficients 
are non-negative, in order to find the maximal overlap we can make the hypothesis that the closest 
separable state is of the form 

(v^io) + yr^ii)) (v^io) + yr^ii))^^-^ 

We further see that in order for the overlap to be maximal, q must be either 1 or 0. For the 
former case, we can further maximize the overlap to get y(2 — y)/2. For the latter case, the 
maximum overlap is -y/y/2, which is less than \/{2 — y) /2 (as < y < 1). Hence, the state 
■s/y\^G) i ^J1 — y \uk) has the entanglement eigenvalue y^(2 — y)/2. 
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As 1 — A^3^(y) is linear in y and — 21og2 Ainax(y) is convex in y, one gets 

X 2 

Esin<PN{x)) = -, £;iog2(pjv(x)) =log2^— ^, (4.10) 

and one of the optimal decompositions is 

1 ^ 

^^'^^^ = 47vIZ Yl \^c.Ak{x)){tpa,p,k{x)\- (4.11) 

k=l a=± )3=u,v 

The above calculations show that for Pn{x), the entanglement depends on the portion x of the 
GHZ in states \tpa,i3,k{x)) {il^a,f3,k{x)\ and it never becomes zero unless there is no GHZ mixture. 

We conjecture that, for > 4, pn{x) has Er{x) = x, with one closest separable mixed state 
being 

|(|0..0)(0..0| + |i..i)(i..i|) + Pk), 

k=l 

which seems plausible as (|0..0)(0..0| + |1..1)(1..1|) is a closest separable mixed state to |*g)- 
4.4 Concluding remarks 

We have presented analytical results on how much entanglement is bounded in two distinct multi- 
partite bound entangled states. The measure we have used to quantify tlicir entanglement is the 
geometric measure of entanglement (GME), whose construction, similiar to the entanglement of 
formation (Ep), is via convex hull. In contrast to GME, Ef has not been explicitly generalized 
to multi-partite states, and hence is still unavailable for these bound entangled states. However, 
under the circumstances discussed previously, the results of fi'sin^ ^^^^ ^ related quantity, 
E'logj, might provide lower bounds on Ep. For the Smolin state, its bound entanglement is as large 
as that of a four-partite GHZ state, whereas that for Diir states is related to the portion of the 
AT-partite GHZ state. For each case, an optimal decomposition is given. Furthermore, we have 
conjectured that the relative entropy of entanglement (E'r) for the Smolin state is unity (proved 
below), whereas we conjecture that -Er for Diir's state is equal to the portion that is A-GHZ. 
For Smolin's state we can establish its Ep, Ejy, E^_ and E^^^i for certain bi-partite partitionings. 
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For example, if we group the four parties ABCD in two, A:BCD, we can write the state as 

= \T,mXi\, (4.12) 



4 

i=0 



with the 3-qubit states of BCD mapped on to the 8-level system (000 0,001 1, Ill 7), 
involving the locally orthogonal and convertible states (by BCD) 

|^o> = (|00) + |17))/^/2, |Xi> = {m + \U))/V2, 
1X2) = (|05) + |12))/V2, 1X3) = (|06) + |11))/V2. 



In order to find the entanglement of this bi-partite state (in C^), we need to consider the 
entanglement of the general (properly normalized) pure state 

that appears in the pure-state decompositions. In fact, regardless of the values of the Xj's, this 
pure state has a reduced density matrix (tracing over BCD) of the form (|0)(0| + |1)(1|) /2. This 
shows that p^^-^CD ^^^^ _ -^^ E^^^2 = 1/2, and -E'iog2 = 1- In feet, there is a general result due 
to Horodecki and co-workers [74] that = E-p for mixture of locally orthogonal bi-partite states, 
e.g., (8) C^™ states that are derived from mixing Bell- like states 



|^^) = (|0,fc)±|l, 2m-A;-l) )/V2, (4.13) 

having distinct k^s, where k = 0,1, . . . ,m — 1. As Ej) < E^ < E-p, we have that Eji{p^''^'~^^) = 1 
as weh. What about the original four-partite state ^^^^07 £'j^(pABCD-j > Eji{p^'-^'^^) , we have 
Epi{p^^'^^) > 1. But we also have that £'r(p^^^^) < 1, as our previous conjecture gives at least 
an upper bound; we thus have that £'r(/9^^^^) = 1 and the conjecture is proved. Naively, we 
expect that any arbitrary p^^^D j^g^g greater entanglement than p^-^CD^ However, for the Smolin 
state, they have the same entanglement as quantified by both GME and the relative entropy of 
entanglement. 

Although Diir's bound entangled state violates a Bell inequality, it has nonzero negativity under 



certain partitionings. One may raise the question: Does there exist a bound entangled state that has 
positive PT (PPT) under all partitionings but that still violates a Bell's inequality? For example, 
does a UPB bound entangled state [65] violate a Bell inequality? We shall see shortly that the 
answer is "No", at least for the three different Bell inequalities [69 \ \71 \ [72] mentioned earlier. Acm 
has shown [75] that if an A^-qubit state violates a two-setting Bell inequality then it is distillable 
under certain bi-partite partitioning. Using the results of Refs. \76\ [77] regarding distillability, we 
can repeat the same analysis for the other two inequalities [Til [72] and indeed obtain the same 
conclusion. We analyze this as follows. 

It was shown by Diir and Cirac [77] that an arbitrary A^-qubit state p can be locally depolarized 
into the form 

2]V-i_i 

p^ = A+|^'o+)(*+|+Aol*o)(*ol+ E A,(|*+)(*+| + |^7)(^-|), 

while preserving = (^'g and Xj = + (^'J l/jl^fT), where j^-Q ) = (jO'^^) ± 

|l®^\)/v/2, and the 



y's are GHZ-like states, i.e., the states in (|4.13p . unfolded into qubit 
notation. Normalization gives the condition 

j 

Now define A = Aj" — Aq , which we assume to be non-negative (w.l.o.g). The condition that there 
is no bi-partite distillability for some bi-partite partitioning Pj is [76) 



2Xj > A. 



Assuming non-distillability for all bi-partite splittings, we have 
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As Aq + Aq > a, we have further that 

1 - A > (2^-1 - 1)A. (4.14) 

For the Mermin-Klyshko-Beh inequahty, violation imphes A > 1/2^^^^^/^. For the three-setting Beh 
inequahty considered in [71j, violation implies A > \/3(2^/3^). For the functional Bell inequality 
in |72j . violation implies A > 2 {2^ /n^). One can easily check that the three Bell inequalities con- 
sidered are inconsistent with non-bipartite-distillability condition, Eq. ()4.14p . Hence, the violating 
of these three Bell inequalities implies the existence of some bi-partite distillability. 

This bi-partite distillability then implies a negative PT (NPT) under that bi-partite partitioning 
according to Horodecki and co-workers [T9]. Hence, violating these Bell inequalities implies NPT 
under certain bi-partite partitioning. Said equivalently, if an A^-qubit state has PPT under all 
bi-partite partitionings then the state never violates these Bell inequalities. This seems to suggest 
that PPT bound entangled states are truly bound in nature that cannot give deviation from local 
theories. 
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Chapter 5 

Global entanglement and quantum 
criticality in spin chains 

5.1 Introduction 

Entanglement has been recognized in the past decade as a useful resource in quantum information 
processing. Very recently, it has emerged as an actor on the nearby stage of quantum many- 
body physics, especially for systems that exhibit quantum phase transitions [THl [79l [801 ED 182] . 
where it can play the role of a diagnostic of quantum correlations. Quantum phase transitions [83] 
are transitions between qualitatively distinct phases of quantum many-body systems, driven by 
quantum fluctuations. In view of the connection between entanglement and quantum correlations, 
one anticipates that entanglement will furnish a dramatic signature of the quantum critical point. 
On the other hand, the more entangled a state is, the more useful it is likely to be as resource 
for quantum information processing. It is thus desirable to study and quantify the degree of 
entanglement near quantum phase transitions. By employing entanglement to diagnose many-body 
quantum states one may obtain fresh insight into the quantum many-body problem. 

To date, progress in quantifying entanglement has taken place primarily in the domain of bi- 
partite systems. Much of the previous work on entanglement in quantum phase transitions has 
been based on bi-partite measures, i.e., focus has been on entanglement either between pairs of 
parties [781 179j or between a part and the remainder of a system [80] . For multi-partite systems, 
however, the complete characterization of entanglement requires the consideration of multi-partite 
entanglement, for which a consensus measure has not yet emerged. 

Singular and scaling behavior of entanglement near quantum critical points was discovered in 
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important work by Osterloh and co-workers [79j . who invoked Wootters' bi-partite concurrence j22j 
as a measure of entanglement. The drawback of concurrence is that it can deal with only two spins 
(each with spin-1/2) even though the system may contain an infinite number of spins. Although 
attempts have been made to generalize concurrence to many spin-1/2 systems via the time reversal 
operation, the generalized concurrence loses its connection to the entanglement of formation |84j . 

Another approach is to consider the von Neumann entropy of a subsystem of L spins with the rest 
N — L spins of the system. It is found that for critical spin chains the entropy scales logarithmically 
with the subsystem size L for N ^ oo, with a prefactor that is related to the central charge of 
the corresponding conformal theory [80\ I81j . However, the entanglement addressed in this case 
is not truly many-body, but only between a subsystem and the rest of the system, although the 
connection to central charge is interesting. 

Quite recently, Barnum and co-workers [82] have developed an entanglement measure, which 
they call generalized entanglement. Instead of using subsystems they use different algebras and 
generalized coherent states to define the entanglement. They have also applied the generalized en- 
tanglement to systems exhibiting quantum phase transitions. Their approach opens a new approach 
to multi-partite entanglement. However, there is no a priori choice of which algebra, amongst all 
possible ones, is the most natural one to use. 

In addtion to spin chains other models that have been studied by using either the von Neumann 
entropy or the concurrence as the entanglement measure include: (i) the super-radiance model, 
in which many two-level atoms interact with a single- mode photon field |85j : and (ii) the one- 
dimensional extended Hubbard model, in which electrons can hop between the nearest neighbors 
and there are Coulomb interactions among electrons on the same site and with nearest-neighbor 
electrons as well [86]. Verstraete and co-workers [57] have recently defined an entanglement length, 
viz., the distance at which two sites can establish a pure-state entanglement at the cost of measuring 
all other sites. They found that this entanglement length is usually greater than the correlation 
length. All these, including the theme of the present chapter, are aimed at approaching many-body 
problems from different, and hopefully fresh, prespectives. 

In the present chapter, we apply the global measure that we have developed in previous chapters, 
based on a geometric picture; it provides a holistic, rather than bi-partite, characterization of 
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the entanglement of quantum many-body systems. Our focus is on one-dimensional spin systems, 
specifically ones that are exactly solvable and exhibit quantum criticality. For these systems, we are 
able to determine the entanglement analytically, and to observe that it varies in a singular manner 
near the quantum critical line. This supports the view that entanglement — the non-factorization of 
wave functions — reflects quantum correlations. Moreover, the boundaries between different phases 
can be detected by the entanglement. 



5.2 Global measure of entanglement 

We quickly review the global measure that we shall use in the present chapter. Consider a general, 
n-partite, normalized pure state: |^') = J2pi--p^ ^PiP2---pnkpi^6p2^ ■ ' ' ^p"'')- If the parties are all spin- 
1/2 then each can be taken to have the basis {||) , ||)}. Our scheme for analyzing the entanglement 
involves considering how well an entangled state can be approximated by some unentangled (nor- 
malized) state (e.g., the state in which every spin points in a definite direction): |<I>) = 
The proximity of |^') to |$) is captured by their overlap; the entanglement of |^) is revealed by 
the maximal overlap [38] 

A„,ax(^') = max|($|^')|; (5.1) 

the larger A^ax is, the less entangled is 1^*). (Note that for a product state, A^ax is unity.) If the 
entangled state consists of two separate entangled pairs of subsystems, A^ax is the product of the 
maximal overlaps of the two. Hence, it makes sense to quantify the entanglement of \^) via the 
following extensive quantity 

i?iog,(^) = -log2ALx(^), (5.2) 

This normalizes to unity the entanglement of EPR-Bell and A^-party GHZ states, as well as giv- 
ing zero for unentangled states. Finite- entanglement is interesting in the context of quantum 
information processing. To characterize the properties of the quantum critical point we use the 
thermodynamic quantity £ defined by 

£ = lim £]y, (5.3a) 

N~*oo 

£N = N-^Eiog,{^), (5.3b) 
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where is the entanglement density^ i.e., the entanglement per particle. 

5.3 Quantum XY spin chains and entanglement 

We consider the family of models governed by the Hamiltonian 

WxY = - X; + ^-J-lw + hai^^ (5.4) 

where r measures the anisotropy between x and y couplings, h is the transverse external field, 
lying along the z-direction, and we impose periodic boundary conditions, namely, a ring geometry. 
At r = we have the isotropic XY limit (also known as the XX model) and at r = 1, the Ising 
limit. All anisotropic XY models (0 < r < 1) belong to the same universality class, i.e., the 
Ising class, whereas the isotropic XX model belongs to a different universality class. XY models 
exhibit three phases (see Fig. 15.10 : oscillatory, ferromagnetic and paramagnetic. In contrast to the 
paramagnetic phase, the first two are ordered phases, with the oscillatory phase being associated 
with a characteristic wavevector, reflecting the modulation of the spin correlation functions (see, 
e.g., Ref. |88j). We shall see that the global entanglement detects the boundaries between these 
phases, and that the universality class dictates the behavior of entanglement near quantum phase 
transitions. 

Before we solve the entanglement of the XY model, we give perturbative analysis of, as an 
illustration of how entanglement arises and vanishes, the Ising model in a transverse field (viz. 
r = 1) 

N 

n = -Y,{afaf^, + ha^). (5.5) 
1=1 

At h = the ground state is that with all spins pointing up in the x-direction f^— > • • • — >) or 
down |— <— • • • <— ) , which is manifestly unentangled. The ground state can be any superposition of 
(]—>—> • • • and |— ^ • • • <— ) when the Z2 symmetry is not spontaneously broken. For example, 
the states (j->— > • • • — >) =b|— <— • • • <— ))/-v/2 are actually the two lowest levels obtained from solving 
the models using Jordan- Wigner and Bogoliubov transformations and they both have -E'logj = 1- 
(For small h the entanglement rises quadratically in the case of unbroken symmetry instead of 
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quartically, as we shall show shortly.) We shall see later that whether or not we use a broken- 
symmetry state has no effect in the thermodynamic limit. For small h (i.e., h <^ l/\fN) one can 
obtain the ground state by treating the haf terms as perturbations. Take the ground state at /i = 
to be ]— >— > • • • — >) . Then first-order perturbation theory for the ground state gives 

h 



J- ^ 4 



(h-----)+2EH----.----))- (5.6) 



Using the method described in Sec. 12.21 we obtain -E'iog2 ^ N{N — l)/i^/32 to leading order in h. At 
h = oo the ground state is a quantum paramagnet with all spins aligning along the external field: 
ft . . . t), and once more is unentangled. To 0{l/h) perturbation theory gives (treating afafj^i 
terms small) 



1 , N V 4/i 
^ + W * 

for which -Eiogg ^ N/{\&h?'), to leading order of 1/h. The quantum phase transition from a 

ferromagnetic to a paramagnetic phase occurs at /i = 1 [83|. The two lowest levels, which we 

denote by |^i/2) ^-iid |^'o) (for reasons to be explained later) are, respectively, the ground and first 

excited states, and they are asymptotically degenerate for < /i < 1 when oo. 

As is well known [83\ [88l 189] . the energy eigenproblem for the XY spin chain can be solved 

via a Jordan- Wigner transformation, through which the spin degrees of freedom are recast as 

fermionic ones, followed by a Bogoliubov transformation, which diagonalizes the resulting quadratic 

Hamiltonian. 

The Jordan- Wigner transformation that we shall make from spins (c's) to fermion particles 
(c's) is 

af = 1-24q, (5.8a) 

i-l 

af = l[{l-2c]c,){c, + 4), (5.8b) 

j-i 

= -^11(1 -24^^) (^^-4)- (5-8c) 

One has to pay attention to the boundary conditions that are to be imposed on the c's. Although 
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periodic in the cr's, one cannot simply take 

N 

= n (1 " ^4 (^^+1 + ^Wi) = ^1 = (ci + 4) : (5-9) 

and conclude that either (i) IljLi (l ~ 2ct Cj) = 1 and cjv+i = ci, or (ii) IljLi (l ~ 2ct Cj) = — 1 and 
cn+1 = — ci, neither of which are correct if one wishes to obtain the correct spectrum and eigenstates 
for arbitrary finite N. Instead one should impose that (as the {N + l)-th site is identified as the 
first site) 

f^AT^AT+l = (5-10) 

which then leads to 

N 

{cN + 4) (cjv+1 + cjv^i) = - n (1 ~ ^4 ^i) (^^ + ^n) (ci + 4) • (5-11) 

i=i 



The two possible conditions that satisfy this are either (I) Ylf=i (l — 2cj Cj) = —1 and cat+i = ci, 
or (II) YljLi (l — 2cJ Cj) = 1 and cn+i = —c\. The operator 



N 

JJ (1 - 24 Cj) = e^" '^^ (5.12) 

counts whether the total number of particles is even (+1) or odd (—1). For c's that are periodic, 

the number is odd, whereas for antipcriodic c's, this number is even. 
To incorporate these two boundary conditions on the c's, we take 



1 

N 



m=0 



where 6 = for periodic c's; b = 1/2 for anti-periodic c's. (This explains why we label the lowest 

two states by 1^*1/2) S'^d l^o)-) The momentum index m ranges from to — 1. In terms of these 
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fermion operators the Hamiltonian becomes 



H = -Nh-Y^ I [2 cos |(m + 6) - 2/.] i^^i^ + ir sin |(m + 6) [ct^ciJL^.,, + ct^tc^L-^J | 

m=0 J 

(5.14) 

Upon using the Bogoliubov transformation 



with 



,an2«) = r,mM|±ii/(ft_c„,M!^), (5.16) 
one arrives at the diagonal the Hamiltonian: 

N-l 

7^ = -AT/. + (5.17a) 

m=0 

41 = 2y'(A-cosM^)^,.,i„.M!!L±i4, (5.1Tb) 

except for the special case that £q°^ = 2{h — 1). 

We remark that we have not left out any constant in diagonalizing the Hamiltonian in either 
case, so the energy spectrum is exact. For each value of b the diagonalization gives 2''^ energy 
eigenvalues, so there are 2^+1 in total. Half of them are spurious. In determining the correct 
2^ states from the 2-^+^ solutions, one has to impose a constraint from the boundary conditions. 
Namely, in case I there can be only odd number of fermions, whereas in case II there can be only 
even number of fermions. 

For 6 = 0, viz, the odd-number-fermion case, the lowest state |^o) is such that {'Jm^jm^) = 
except that (70°'* ^7o'^'*) = 1- Its energy eigenvalue is 



Af-l 



4°' (., ft) = (fc - 1) - ^ (ft - cos % sin' (5.18) 

m=l 

For b = 1/2, namely, the even-number-fermion case, the lowest state |*i/2) is such that 



96 



(7m''^^^7m'^^^) = for all m. Its eigen-energy is 

Er\r, k) = - g y(.-co.M^)V..,i,.M!!!±iZ?). (5.19) 

m=0 

We see that, as ^ oo, the above two energy levels are degenerate for h < 1. Furthermore, as 
N ^ oo the difference between the two energy levels becomes 

E^°\r,h) - Ej^^^^\r,h) = 2{h - l)e{h - 1), (5.20) 

where Q{x) = 1 if x > and zero otherwise. The way the energy gap vanishes as /i ^ 1^ gives a 
relation betweem two exponents 

zi^ = l; (5.21) 

z is the dynamical exponent (defined via the vanishing of energy gap A ~ \h — hcY^) and u is the 
correlation-length exponent (defined via Lc ~ |/i — hc\~^). 

Having found the lowest two eigenstates, the quantity Amax of Eq. (j5.ip — and hence the entanglement — 
can be found, at least in principle. To do this, we parametrize the separable states via 

\^)=® [cos(e./2)t), + e^<^'sin(e./2)i),], (5.22) 

where |t/J,) denote spin states parallel/antiparallel to the z-axis. Instead of maximizing the overlap 
with respect to the 2N real parameters {Ci, (/>«}, for the lowest two states it is adequate to appeal 
to the translational symmetry of and the reality of the ground-state wavefunctions. Thus taking 
= i and 0i = we make the Ansatz: 

|^>(e)) = e~^2^^-i"^' ITT . . . T) (5.23) 

for searching for the maximal the overlap Ajnax(^) [SD]. This form shows that this separable state 
can be constructed as a global rotation of the ground state at /i = oo, viz., the separable state 
ITT • • • T)- III this particular limit the boundary condition on the c's is irrelevant, as the dominant 
term in the Hamiltonian is — ^ • h{l — 2CjCj). 
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The energy eigenstates are readily expressed in terms of the Jordan- Wigner fermion operators, 
and so too is the family of the Ansatz states |#(^)). By working in this fermion basis we are able 
to evaluate the overlaps between the two lowest states and the Ansatz states. With |^'o) (|^i/2)) 
denoting the lowest state in the odd (even) fermion-number sector, we arrive at the overlaps 

{^,{r,h)\m)=f^N\0 n [^oseg\r,h)cos\cm+smei^\r,h)sm\C/2)cot{kl^]^/2)] , 

m=l-2b 

(5.24) 

with 

€]N^'^im + b), te.n2e^^\r,h) ^ rsmk^^]j,/ih-cosk^^]j,); (5.25a) 
fW^HO ^ 1, fS\0 ^ cos(e/2), (N even); (5.25b) 

fS^^HO ^ cos(e/2), fS\c) ^ VNsmiC/2), {N odd); (5.25c) 

where 6 = 0,1/2 and m E [0, — 1] is the (integer) momentum index. The above results are exact 
for arbitrary A^, obtained with periodic boundary conditions on spins rather than in the so-called 
c-cyclic approximation [89j. Given these overlaps, we can readily obtain the entanglement of the 
ground state, the first excited state, and any linear superposition, cos a \ ^o) + sin a \ ^i) of the two 
lowest states, for arbitrary (r, h) and A^, by maximizing the magnitude of the overlap with respect 
to the single, real parameter ^. For the derivation of the above results, see Appendix [Gl 

The formulas [in Eqs. (I5.24P and (j5.25p ] contain all the results that we shall explore shortly. By 
analyzing the structure of Eq. ()5.24p , we find that the global entanglement does provide information 
on the phase structure and critical properties of the quantum spin chains. Two of our key results, 
as captured in Figs. 15.11 and 15. 2^ are: (i) although the entanglement itself is, generically, not 
maximized at the quantum critical line in the (r, h) plane, the field- derivative of the entanglement 
diverges as the critical line h = 1 is approached; and (ii) the entanglement vanishes at the disorder 
line + /i^ = 1, which separates the oscillatory and ferromagnetic phases. 

As is to be expected, at finite A^ the two lowest states |^'o) and |*I'i/2) featuring in Eq. ()5.24p 
do not spontaneously break the Z2 symmetry. However, in the thermodynamic limit they are 
degenerate for h < 1, and linear combinations are also ground states. The question then arises 
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Figure 5.1: Entanglement density (upper) and phase diagram (lower) vs. {r,h) for the XY model 
with = 10^ spins, which is essentially in the thermodynamic limit. There are three phases: O: 
ordered oscillatory, for + /i^ < 1 and r 7^ 0; F: ordered ferromagnetic, between + /i^ > 1 and 
/i < 1; P: paramagnetic, for h > 1. As is apparent, there is a sharp rise in the entanglement across 
the line h = 1, which signifies a quantum phase transition. The arc /i^ + = 1, along which the 
entanglement density is zero (see also Fig. 15. 2p . separates phases O and F. Along r = lies the 
XX model, which belongs to a different universality class from the anisotropic XY model. 
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Figure 5.2: Entanglement density and its /i-derivative (inset) for the ground state of three systems 
at = oo. Sohd hne: Ising (r = 1) limit; dashed line: anisotropic (r = 1/2) XY model; dash- 
dotted line: (r = 0) XX model. For the sake of clarity, the XY-case curves are shifted to the 
right by 0.5, indicated by the arrow. For the r = 1/2 case, at the entanglement density vanishes 
at /i = -y/l — r^, which is a general property for the anisotropic XY model. Note that whilst the 
entanglement itself has a nonsingular maximum at /i ~ 1.1 (Ising), h ~ 1.04 (r = 1/2 XY), and 
h = Q (XX), respectively, it has a singularity at the quantum critical point at /i = 1, as revealed 
by the divergence of its derivative. 

as to whether linear combinations that explicitly break Z2 symmetry, i.e., the physically relevent 
states with finite spontaneous magnetization, show the same entanglement properties. In fact, we 
see from Eq. (j5.24p that, in the thermodynamic limit, overlaps for |^o) ai^d I ^1/2) ^'^^ identical, up 
to the prefactors and f^^'^^ ■ These prefactors do not contribute to the entanglement density, 
and the entanglement density is therefore the same for both |^'o) and |^'i/2)- It further follows 
that, in the thermodynamic limit, the results for the entanglement density are insensitive to the 
replacement of a symmetric ground state by a broken-symmetry one. 



5.4 Entanglement and quantum criticality 

Before we discuss the thermodynamic limit of the entanglement density, we compare the entan- 
glement obtained via the results in Eq. (15.241) and that obtained via numerically diagonalizing the 
Hamiltonian and calculating the maximal overlap. In Figure 15.31 the results via each method are 
shown for the Ising case (r = 1) with small numbers of spins {N = 13 through 22), it is seen that 
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Figure 5.3: Entanglement vs. magnetic field h for ground (upper panel) and first excited (lower 
panel) states with small N = 13, 16, 19, 22 (from bottom to top) for the transverse Ising model 
(r = 1). The numerical results are shown as discrete points whereas the analytical results are 
shown as lines. This demonstrates that the analytical results are exact, even for small N, both 
even and odd. 

our analytical results are exact even for small A^, both even and odd. 

From Eq. (|5.24p it follows that the thermodynamic limit of the entanglement density is given 

by 

£{r,h) = max / d/i In [cos 0(/i, r, /i) cos^(^/2) + sin6'(^, r, /i) sin^(,^/2) cot vr^l , (5.26) 

In 2 € Jo L 

where tan 20(//, r,h) = r sin 2TTfi/{h — cos 27r^). 

Figure 15.21 shows the thermodynamic limit of the entanglement density £{r,h) and its h- 
derivative in the ground state, as a function of h for three values of r, i.e., three slices through the 
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surface shown in Fig. 15.11 As the r = 1 shce shows, in the Ising hmit the entanglement density is 
smah for both small and large h. It increases with h from zero, monotonically, albeit very slowly 
for small /i, then swiftly rising to a maximum at /i ~ 1.13 before decreasing monotonically upon 
further increase of /i, asymptotically to zero. The entanglement maximum does not occur at the 
quantum critical point. However, the derivative of the entanglement with respect to h does diverge 
at the critical point /i = 1, as shown in the inset. The slice at r = 1/2 (for clarity, shifted half a 
unit to the right) shows qualitatively similar behavior, except that it is finite (although small) at 
/i = 0, and starts out by decreasing to a shallow minimum of zero at /i = \/l — r^. By constrast, 
the slice at r = (XX) starts out at /i = at a maximum value of 1 — 27c/(7r ln2) ~ 0.159. (where 
7c ~ 0.9160 is the Catalan constant), the globally maximal value of the entanglement over the 
entire (r, h) plane. For larger h it falls monotonically until it vanishes at /i = 1, remaining zero for 
larger h. 

We find that along the line r'^ + h? = 1 the entanglement density vanishes in the thermodynamic 
limit. In fact, this line exactly corresponds to the boundary separating the oscillatory and ferromag- 
netic phases; the boundary can be characterized by a set of ground states with total entanglement 
of order unity, and thus of zero entanglement density. The entanglement density is also able to 
track the phase boundary [h = 1) between the ordered and disordered phases. Associated with 
the quantum fluctuations accompanying the transition, the entanglement density shows a drastic 
variation across the boundary and the field-derivative diverges all along h = \. The two boundaries 
separating the three phases coalesce at (r, /i) = (0,1), i.e., the XX critical point. Figures ET] and 
15.21 reveal all these features. 

In Appendix IH.ll we analyze the singular behavior of the field-derivative of the entanglement 
density ()5.26p in the vicinity of the quantum critical line, and we find that the asymptotic behavior 
(for r 7^ 0) 

d£ 1 

From the arbitrary- results (I5.24P of the entanglement we analyze the approach to the thermo- 
dynamic limit, in order to develop further connections with quantum criticality. We focus on the 
exponent which governs the divergence at criticality of the correlation length: Lc ~ \h — 1|~''. 
To do this, we compare the divergence of the slope dS^ /dh (i) near h = 1 {a.t N = oo), given 
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Figure 5.4: Finte-size scaling. dS^ /dh\h^^^ j^ vs. IniV, for N = 10^^,3 x 10"^, 5 x 10'*, 8 x 10^, and 
10^ (points) with r = 0.1. The solid line represents the fit dS^ /dh\h^^^^ k. 2.30 In — 6.95. 

above, and (ii) for large N at the value of h for which the slope is maximal (viz. /imax.Af)) i-e., 
dE^ / dh\h^^^ ~ 0.230r^-'^ InA^ + const., obtained by analyzing Eq. (j5.24p for various values of r; 
see Fig. 15.41 for the example of r = 0.1 case. Then, noting that (27rln2)^* ~ 0.2296 and that the 
logarithmic scaling hypothesis [9T] (or see Appendix[Il) identifies v with the ratio of the amplitudes 
of these divergences, 0.2296/0.230 ~ 1, we recover the known result that v = 1. Moreover, from 
Eq. (j5.2ip we can extract the value of the dynamical exponent: z = 1. 

We derive in Appendix IH.2I the divergence behavior of the field-derivative of the entanglement 
density for the isotropic (r = 0) case. Compared with r ^ case, the nature of the divergence of 
d£/dh at r = belongs to a different universality class: 

(0, h) . - ih^l). (5.28) 

From this divergence, the scaling hypothesis, and the assumption that the entanglement density is 
intensive, we can infer the known result [83] that the critical exponent u = 1/2 for the XX model. 
Moreover, from Eq. (j5.2ip we can extract the value of the dynamical exponent z = 2 for the XX 
model. 

In keeping with the critical features of the XY-model phase diagram, for any small but nonzero 
value of the anisotropy, the critical divergence of the entanglement derivative is governed by Ising- 
type behavior. It is only at the r = point that the critical behavior of the entanglement is 
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governed by the XX universality class. For small r, XX behavior ultimately crosses over to Ising 
behavior. 

We have mentioned that along the disorder line, r^ + /i^ = 1, the entanglement density vanishes. 
One extreme limit is the Ising case, i.e., r = 1 and h = 0, where the ground state is either ]— >— > • • • — >) 
or |— <— • • • <— ), both of these being unentangled. Any superposition of them is also a valid ground 
state, but it has entanglement of order unity. In the thermodynamic limit the entanglement per 
spin is identically zero. Is this a general feature along the disorder line? Before we establish this, 
recall that the energies of the lowest two levels are given in Eqs. (I5.18|) and ()5.19p . Evaluating them 
at = 1 — /i^, we immediately find that both are — A^. 

Now let us evaluate the expectation value of the Hamiltonian (15. 4p with respect to a separable 
state with all spins pointing in the same direction: 



Denoting x = {(7^), y = (c^), z = (cr^), we find that the above expression achieves its minimum 
value — A at + /i^ = 1 when 



Therefore, the separable state satisfying the above conditions is the ground state, and there is a Z2 
degeneracy. Hence, along the disorder line the entanglement density vanishes. 

5.5 Concluding remarks 

In summary, we have quantified the global entanglement of the quantum XY spin chain. This 
model exhibits a rich phase structure, the qualitative features of which are reflected by this entan- 
glement measure. Perhaps the most interesting aspect is the divergence in the field-derivative of 
the entanglement as the critical line {h = 1) is crossed. The behavior of the divergence is dictated 
by the universality class of the model. Furthermore, in the thermodynamic limit, the entanglement 
density vanishes on the disorder line (r^ + h? = 1). The structure of the entanglement surface, as 




(5.29) 




(5.30) 
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a function of the parameters of the model (the magnetic field h and the coupling anistotropy r), is 
surprisingly rich. 

We close by pointing towards a deeper connection between the global measure of entanglement 
and the correlations among quantum fluctuations. The maximal overlap (jS.ip can be decomposed 
in terms of correlation functions (see Sec. 12.2.^ : 

■y Ji^ ^ TV 

^max = ^ + ^ {^^' "^^^ + 2 '^1 ® '^i) + • • • }' 

J— 2 

where translational invariance is assumed and the Cartesian coordinates of r can be taken to be 
(sin ^, 0, cos ^) . The two-point correlations appearing in the decomposition are related to a bi-partite 
measure of entanglement, namely, the concurrence, which shows similar singular behavior [TS] to 
Eq. (|5.27p . It would be interesting to establish the connection between the global entanglement and 
correlations more precisely, e.g., by identifying which correlators are responsible for the singular 
behavior in the entanglement and how they relate to the better known critical properties. 
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Appendix A 



Local hidden variable theories and 



Bell-CHSH inequality 



Local hidden variable (LHV) theories were conceived as alternatives for quantum mechanics and 
were constructed in order to explain the predictions of quantum mechanics in a statistical way 
without invoking the nonlocal features that quantum mechanics allows. The most natural quantity 
for two observers at a distance to measure is the correlation. Local hidden variable theories dictate 
that the outcome of some measurement is predetermined by the combination of measurement setting 
and some unknown local hidden variable, and that the result at one site should not be influenced 
by that at the other. Suppose A and B are operators to be measured at the two different sites, 
respectively. The correlation predicted by the LHV theories is 



where A{a, A) = ±1 and B{b, A) = ±1 are predetermined results for the measurement settings a for 
A and b for B depending on the local hidden variable A; p{X) is the distribution for the local hidden 
variable. Locality requires that the outcome A{a, A) does not depend on b and that of B{b, A) does 
not depend on a. For two different settings at each site one arrives at the inequality 




(A.1) 



\EL{a, b) + EL{a, b') + EL{a', b) - EL{a', b')\ < 2. 



(A.2) 
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We shall see shortly that quantum mechanics can violate this inequality. To be specific, the oper- 
ators to be measured are the Pauli operators a. The correlation 



EQ{a,h) = {a -0,(^5 -h), (A.3) 

where a and b are unit vectors indicating the directions of the Stern-Gerlach apparata at A and 5, 
respectively. Define 

2B = a- a®a-h + a- d®a-U + a-a'®a-h-a-a'®a-b'. (A.4) 
It is straightforward to see that for a separable state ps = Yli PiPa Pb 

max \2Tt{Bps)\ < 2, (A.5) 

a,a' ,b,b' 

namely, it never violates the inequality. But for a singlet state (|t|) — |||))/\/2, we have 

max 1(25)1 = 2V2. (A.6) 

a,a',b,b' 

This can be achieved for the settings 6a = 0, 6*^ = it/2, 6), = 7r/4, and 6'^^ = — vr/4, where the angles 
are measured from the z-axis in the z — x plane. 

Gisin [4j was the first to show that for any entangled pure state, e.g., a |TT) + /? with a(3 / 0, 
the inequality is also violated. Horodecki and co-workers [92| derived the maximal violation for any 
two-qubit mixed states. It was then clear that there exist many mixed states that are entangled 
but do not violate the CHSH inequality. The question remains open whether there exists a Bell-like 
inequality that is necessary and sufficient for a state being entangled. 
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Appendix B 



Schumacher's quantum data 



compression 



B.l Quantum Data Compression 

To communicate quantum information by directly transmitting qubits may be costly. The idea of 
quantum data compression (QDC) is to ask the question whether we can compress the message 
into fewer qubits so as to minimize the cost of transmission. Schmacher has provided an answer to 
achieve this |21j . One excellent review of quantum data compression is by Preskill [9^, which we 
follow here. 

Let us go straight to the procedure of QDC. Suppose Alice needs to communicate with Bob 
through some noiseless quantum channel as efficiently as possible, that is, she hopes to compress 
her message using as few qubits as possible. The message consists of letters represented by some 
states \4>x)- Since on average, the frequency of each letter's appearance may not be equal but is 
some probability px-, the message can be said to be drawn from an ensemble of states: 



so each letter has a density matrix p = ^xP^\^^){^^\- ^'^^ ™- following, the 

lowest number of qubits per letter needed to encode is set by the von Neumann entropy S{p) = 
— tr(/9log2 p). If we try to compress into fewer qubits, the fidelity of compression will be ruined. 
If the total length of the message is n, then the message has a density matrix which is a direct 
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product of n letter density matrices 



p®" = p O ■ ■ ■ ® p . (B.2) 

n p's 

The procedure for quantum data compression goes as follows, 

1) Diagonalize p. Work in the orthonormal basis in which p is diagonal. If p has eigenvalues 
(arranged decreasingly) Ai > A2 > • • • > Ad(the number d depends on whether the state is a 
qubit, tri-bit, or d-bit), then p®" has eigenvalues of the form (i.e., the eigenvalues are obtained by 
choosing n values from Ai, A2, . . . , Xa) 

A(W) = ^A,^^=A^A2^^•••A^^ (B.3) 

i=l 

where Yli — ^^'^ each eigenvalue X{{ki}) occurs N{{ki}) = n\/{ki\k2l • • • /c^!) times. We will 
restrict ourselves to qubits, i.e., d = 2: \{ki,k2 = n — ki) = A^^A2^, and N{ki,k2) = C^^. 

2) Given a set of tolerances S (tolerance for using slightly more qubits than the asymptotically 
optimal case) and e (tolerance for not projecting onto the "typical" subspace), find the typical 
subspace A and its dual subspace A-*-: 

2a) First find out the smallest number D(n) of necessary largest eigenvalues (suppose they are 
An,i > K,2 > ■ ■ ■ > K,d) and corresponding eigenvectors(|A„,i), |A„,2), • • • AK,d)) of p*^" such 
that the sum of these eigenvalues is larger than 1 — e, but the value D{n) may be larger than 
2"'('S'(p)+<5) . Increase n and repeat the step until when n > uq 

£>(«)< 2"(^(^)+^), (B.4) 

where no is the smallest number such that the above inequality is satisfied. Note that it is sufficient 
to use at most some n{S{p) + S) qubits to represent any state in A. This means there exists a 
unitary transformation U which takes any state in A to 

U\(i)A) = I '/'compressed) I Orest)' (B-5) 
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where l^compressed) ^ state of n{S{p)+6) qubits, and |Opgg|;) is a state \0)(d- • • |0) of {n—n{S+6)) 
qubits. 

2b) Those eigenvectors corresponding to the first D{n) largest eigenvalues span a typical sub- 
space A, the remaining spanning a dual subspace A-*-. Note this division into two subspaces can be 
represented by a projection operator E which projects onto A and the complement of which 1 — E 
to A"*-. The condition that the sum of eigenvalues of eigenvectors in A is larger than 1 — e can be 
rewritten as 

tr(p"£;) > 1 - e. (B.6) 

This means states in A have much higher overlap with any state drawn from the ensemble than 
those in A-*-. 

3) Prepare the input state \ip) = |0i) • • • |</>n), where \<j)i) belongs to the ensemble in the Eq. (lB.l|l . 

Make the unitary transformation U on and measure the state of the last (n — n{S + 6)) qubits 

mentioned above. If the result is |Opgg|^), Alice successfully compresses |^) onto IV'compressed^l^rest)' 

and she simply sends IV'compressed) *° Bob. On the other hand, if Alice gets the results other than 

lOrpot), she fails to compress her message and the best she can do is send a state |0' i) 
I lesi/' f to I compressed' 

which is the compressed state corresponding to the largest eigenvector |Ari,i) in A, 

= |0;ompressed)|Orest)- (B.7) 

We note that the input state {ip) has much higher overlap with states in A than any other states 
in A-*-, the result for Alice to get |Opgg|;) is of high probability (larger than 1 — e). 

4) Bob, after receiving iV'compressed)' ^-PP^nds |Opgg|;) to it, and applies the inverse unitary trans- 
formation U~^. On average. Bob receives a density matrix 

= E\i,){i;\E + |A„,i)(A„,iKV^|(l - Em. (B.8) 

The averaged fidelity F of this procedure over the ensemble of possible messages {\i^i),p'i} can 
be shown to be larger than 1 — 2e. It can also be shown that if we try to compress the message 
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into n{S{p) — S) qubits, the fidelity will be arbitrarily small for sufficient large n. 

Finally, we note quantum data compression cannot compress messages drawn from a com- 



B.2 An Example 

The example we will discuss shortly is for a small n. Prom previous discussion, we know for any 
given 6 and e, we can always find a number no, such that for any n > uq, the procedure succeeds 
with the prescibed tolerances. In fact, for a given n, as e becomes smaller, the necessary S increases, 
which means we need more qubits to compress, as can be seen from the average number of qubits 
necessary to encode is n{S{p) + 6). On the other hand, for a given n, as 6 decreases (we require 
fewer qubits to encode), e increases, which means the average fidelity decreases, as can be seen 
from F > 1 — 2e. Hence, there is some tradeoff between 6, and e for a fixed finite n. 

Suppose the ensemble consists of {(\H),pH = 5),(|-D),]3d = 5)}, where \H) is the state of 
horizontal polarization while \D) is 45°, 



pletely (maximally) mixed state, since S{p completely mixed) ~ 




(B.9) 



The density matrix is 




(B.IO) 



where the matrix is in \H) and \V) basis, and has two eigenvectors and eigenvalues 




(B.ll) 
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The von Neumann entropy S{p) is 



S{p) = -Xq log2 Aq - Aq log2 Xq ^ 0.60088, 



(B.12) 



so the minimal number of qubits per letter needed to encode is about 0.6009. 

Suppose Alice needs to send 3 letters to Bob, but she can afford only two qubits. Since 3*S{p) « 
1.8, it's possible to compress 3 letters using only 2 qubits with high fidelity. Note this means S 0.2 
and D{n = 3) = 2. The eigenvalues and eigenvectors of are 

Ai = cos^ I , A2 = A3 = A4 = cos^ I sin I , 
A5 = Ae = A7 = cos f sin^ As = sin^ f . 



|1) = \QQQ), |2) = \QQQ), |3) = \QQQ), |4) = \QQQ), 



|5) = \QQQ), |6) = \QQQ), |7) = \QQQ), |8) = \QQQ). 



(B.13) 



The subspace A is spanned by {|1), |2), |3), |4)}, while its dual subspace A"*- is by {|5), |6), |7), |8)}, 
and 



Pa = tT{p®^E) = ^Xif^ 0.9419 

i=l 

8 

= tr(p®3(i _ E)) = ^X^^ 0.0581. (B.14) 
Alice and Bob both agree on the form of the unitary transformation that they will use, 



U 







\HHH) ^ 








' \HHV) ^ 


|2) 




\HVH) 


U 


|6) 




\HVV) 


|3) 




\VHH) 


|7) 




\VHV) 


Vi4) J 




\ \VVH) , 




I 18) J 




\ \vvv) J 



(B.15) 



this transformation is unitary since it is simply the transformation between two sets of orthonormal 
bases of 3 qubits. 
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Alice prepares her message in a state {ip), which can be expanded in the basis of |1) • • • |8), 

8 

|^) = ^a,|i), (B.16) 

i=l 

where from Eq. (lB.14]l we have 

4 8 

^|ai|2 = PA» J^|ai|' = ^A- (B.17) 



i=l i=5 



Then Ahce apphes the unitary transformation U on followed by a measurement on the third 
qubit. If the result is \H), she sucessfully projects \ip) into the likely subspace A. At this stage, the 
total state is 

ai\HHH) + a2\HVH) + a^\VHH) + ai\VVH) = iV'compressed)!^)' (^-18) 

where IV'compressed) ^ ai\HH) + a2\HV) + a^\VH) + a4|yy). She simply sends this two-qubit 
state IV'compressed) *° Bob. Upon receiving iV'compressed)' -^^^ appends a third qubit \H) to it, 
and does the inverse transformation to get 

4 

W) = tZ-^dV'compressed)!^)) = E«^l^)' (^.19) 

4 = 1 

which has high resemblance to the initial \tp), that is 

Fi = KV'IV'')!^ = ~ 0-9419. (B.20) 

On the other hand, if, when Alice measures the third qubit and gets \V), she fails to project \ip) 
into A but A-*- instead, the best she can do is send a qubit state \HH). After Bob receives it and 
decompresses it, he gets 

W') = U-WHH)\H)) = \l), (B.21) 
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which has overlap with the initial 



F2 = K^IV'")!^ = kiP = Ai fii 0.6219. 



(B.22) 



The fidelity of this procedure is F = PaFi + P^F2 = 0.9234. 

How good is this? Let us compare it to the case when Alice sends the first two letters without 
compressing and asks Bob to guess the third letter. Since both \H) and \D) from the ensemble 
have higher overlap with \Q) than with \Q), the best guess he can make is \Q). The fidelity of this 
procedure is 



F=^-\{H\Q)\^ + ^\{D\Q)\^ = 0.8535 



(B.23) 



which is smaller than the case when we do compression. 
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Appendix C 

Wootters' formula 



The entanglement of formation defined in Eq. (jl.30p is, in general, difficult to calculate. However, 
for two-qubit systems, Wootters [22j has provided and proved the following formula: 

Ej.{p)=h(^^[i+^i-c{pr]y (c.i) 

where h{x) = —x log2 x — {1 — x) log2(l — x), and C{p), the concurrence of the state p, is defined as 

Cip) = max{0,y^- y^- y^4}, (C.2) 

in which Ai, . . . , A4 are the eigenvalues of the matrix p{cry ® ay)p*[ay ® ay) in nonincreasing order 
and Uy is a Pauli spin matrix. Ep{p), C{p) and the tangle t{p) = C{p)'^ are equivalent measures 
of entanglement, inasmuch as they are monotonic functions of one another. For pure state a |00) + 
6|01) + c|10) + the concurrence C is 2\ad — hc\. 
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Appendix D 

Proof of entanglement monotone 



In this appendix we prove in detail that the geometric measure of entanglement satisfies the criteria 
listed in Sec. 12.31 and hence it is an entanglement monotone. For convenience we list those criteria 
as follows 

CI. (a) E{p)>0; (b) E{p)=0 if /? is not entangled. 
C2. Local unitary transformations do not change E. 

C3. Local operations and classical communication (LOCC) (as well as post-selection) do not 
increase the expectation value of E. 

C4. Entanglement is convex under the discarding of information, i.e., Y2iPi ^(Pi) — ^iYliPi Pi)- 



it is evident that CI and C2 are satisfied, provided that -Epure satisfies them, as it does for -Epure being 
any function of Amax consistent with CI. It is straightforward to check that C4 holds, by the convex 
hull construction. First, consider the case is which p = From the definition (|2.57p of 

E{p), which is the minimum over all decompositions, we have that E{p) < impure (I ^j))- Hence 
we have that E{Y^^pi\tpi){ipi\) < PiE{\ipi) {ipi\) , i.e., C4 is obeyed whenever the deomposition is 
into pure states. Second, allow pi to be mixed. To deal with this case, express pi as its optimal 
decomposition: pi = Y.k9ik\ipik){ipik\, for which E{pi) = Y^,,qikEpuTe{\tpik))- Inserting the above 



From the definition (j2.57p 



E{p) = {coEpnrc){p) = min 
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expression for pi and E{pi) into the left hand side of the sought criterion, and using the pure- 
state result just proved, we find YjiPi^iPi) = J2ikPi 1ikE{\'^ik){'4'ik\) > E{Y,ikPim\'>Pik){ipik\) = 
E{J2iPiPi)- Thus we see that C4 is indeed obeyed. 

The consideration of C3 seems to be more delicate. The reason is that our analysis of whether 
„. not it holds depends on the expheit totn, of E,„.. Fot C3 to hold, ,t is snfficent to show 
that the average entanglement is non-increasing under any trace-preserving, unilocal operation q 
p — > Ylk ^kpyli where the Kraus operator has the form 1 • • • 1 (gi V^^^ ® 1 • • • ® 1 and obeys 
Efc Vk^k = 1. Further more, it suffices to show that C3 holds for the case of a pure initial state, 
i.e., p = 

We now prove that C3 holds for the particular (and by no means un-natural) choice impure = 
-Egj„2. To be precise, for any quantum operation on a pure initial state, i.e., 

we aim to show that 

Y^PkE^n- {yk\i^)/Vp~k) < i^sin^dV')), (D.2) 
k 

where pk = Tr Vfc|'i/')('0l^fc = (^I^J^fclV')) regardless of whether the operation {14} is state-to- 
state or state-to-ensemble. Let us denote by A and the respective entanglement eigenvalues 
corresponding to \ip) and the (normalized) pure state yk\''P) / ^/Pk ■ Then our task is to show that 
"^kPk^k — of which the left hand side is, by the definition of A^, equivalent to 

Vpfe max imiVkm/^f = V max |K4|14|V)f • (D.3) 

k k 

Without loss of generality, we may assume that it is the first party who performs the operation. 
Recall that the condition ()2.9p for the closest separable state 

10) = |a)i® |7)2...„ (D.4) 



^What we mean by unilocal is that the operation is performed by only one of the parties. All general multi-party 
operations can be regarded as a sequence of unilocal operations 
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can be recast as 

2-n(7|V'>l-n = A|a)l. 

Then, by making the specific choice 



(D.5) 



{^k\ = i{a\V^'" /VQ~k) ^ {% (D.6) 
where qk = {oi\V^^^'' V^^^ \a) , we have the sought result 

Y^PkAl = 5^ max \mmm' > a' 5]((a|T^f V«|a)/V¥)' = (D.7) 

k k k 

Hence, the form 1 — A^, when generalized to mixed states, is an entanglement monotone. We note 
that a different approach to establishing this result has been used by Barnum and Linden |37j . 
Moreover, using the result that X^^PfcA^ > A^, one can further show that for any convex increasing 
function fdx) with x G [0, 1], 

J2Pkml)>fciA'). (D.8) 

k 

Therefore, the quantity const. — fc{^) (where the const, is to ensure the whole expression is 
non- negative) , when extended to mixed states, is also an entanglement monotone, hence a good 
entanglement measure. 
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Appendix E 

Three A^-qubit Bell inequalities 



E.l Mermin-Klyshko-Bell inequality 

In Appendix |X] we have discussed the Bell-CHSH inequahty for two qubits. Here we will discuss 
its generalization to N qubits. The discussion here follows Diir [69j. First let us define a Pauli 
operator at arbitrary direction 

aa^=a ■ Ofc, (E.l) 

where is a unit vector. The Mermin-Klyshko-Bell inequality is conveniently defined via Bell 
operators in a recursive way: 




(E.2) 



with Bi = (Tai , B't^ = a^i^ , and i?^ is obtained from by exchanging all with a'^ and vice versa. 
For example, 

B2 = {cTai O-aa + C^a'^ ® ^^aa + C^ai ® C^a'^ " (^a[ ® Cra^)/2 (E.3a) 

It is easy to see that for real x, x', y, y' in the range [—1, 1], we have that \xy+xy' ^x' y—x' y'\ < 2. 
As (Bi) and (i?^) are evidently in the range [—1, 1], local hidden variable theories will then predict 
that |(-B2)| ^ 1; and by induction \{Bk)\ < 1. However, quantum mechanics can violate this 
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inequality. In particular, for the A^-partite GHZ state 



|GHZ) = (|0...0) + |1...1))/V2, 



(E.4) 



I (-Bat) I can achieve the value 2^^ , 

E.2 Three-setting Bell inequality 

There are actually at least two approaches to define a three-setting Bell inequality. These all involve 
selecting three settings for measurement at each site. The first approach is to consider the linear 



possible settings. The goal is to construct an inequality that is satisfied by local hidden variable 
theores whereas it can be violated by quantum mechanics. However, the number of such inequalities 
grows with the number of parties. It is only possible to exhaust all inequivalent inequalities for 
small number of parties, as investigated by Collins and Gisin [94] very recently. Here we shall focus 
on a second approach 3. 

The basic idea of the second approach is that, for two vectors h and if K/ilg)] < H^/lP, this 
means that h ^ q. Specifically, h represents prediction from local hidden variables whereas q 
represents that from quantum mechanics. To be more precise, quantum mechanics predicts that 
for a state pN the average outcome of the observables {0,, i = 1, . . . , N} (with possible outcomes 
being ±1) at locations is 



where is the measurement setting at the fc-th location. A local hidden variable theory predicts 



k=l 

where p{\) is the probability distribution for the local hidden variable A, and IkHk^^) is the 
outcome ascribed by the local hidden variable A for the observable Ok{ik) measured with the 
^It would be interesting to compare the difference between the two approaches. 



combination of the joint probabilities or correlations £^(^1 



. . . , ^n), with being chosen from three 



EQ{ii,...,iN)=TT{Oi---ONPN), 



(E.5) 




(E.6) 
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apparatus setting ^k- 

The inner product between £'(^1, . . . ,£,n) and £"(^1, . . . ,^Ar) is defined as 



{E\E')^ ^ EiCi,...,CN)E'iCi,...,CN), (E.7) 

where each ^ has three different values. If it can be shown that \{Eq\El)\'^ < \\Eq\\'^ then the two 
theories have different predictions. Zukowski and Kaszlikowski |95j considered the following basis 
to be measured: 

|±,e,),^(|0),±|l)J/V2 (E.8) 

for party i. For the N-partite GHZ state they chose three values for ^1: (vr/G, 7r/2, Svr/G), and for 
^f. (0, 7r/3, 27r/3) for z = 2, . . . , A^. They showed that 

\{El\Eq)\ < 2^-1^ < WEgf = 3^/2. (E.9) 

Thus, the two theories, LHV and quantum mechanics, have different predictions. 

E.3 Functional Bell inequality 

The idea of functional Bell inequality [72] is to consider a continuous range of measurement settings, 
instead of a finite number. The inner product between -£(^1, . . . ,S,n) and E'{^i, . . . ,^Ar) for the 
continuous setting is defined as 

N 

{E\E')^ ['[ld^,]E{^i,...,^N)E'{^i,...,CN)- (E.IO) 

i=i 

Similarly, if it can be shown that KEqIEl)]"^ < \\Eq\\'^ then the two theories have different predic- 
tions. Sen and co-workers considered the same measurement basis as in Eq. (jE.Sp . with G [0, 27r]. 
They showed that 

\{El\Eq)\<a'' <\\EQf = {27T)^/2. (E.ll) 

Thus again, the two theories, LHV and quantum mechanics, have different predictions. 

To compare the three-setting and the functional inequalities, we mention that for the Diir's 
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bound entangled state ()4.6p discussed in Sec. 14.31 violation is achieved using the three-setting one 
only when N > 7 whereas for the functional one > 6, the latter being stronger. 
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Appendix F 

Unextendible product bases and 
bound entangled states 

To illustrate the idea of the unextendible product basis and its connection to bound entanglement, 
we shall give an explicit example. The discussion here follows Ref. |65j . product basis is a set S 
of pure orthogonal product states, which span a subspace Tis of "H, a bi- or multi-partite quantum 
system. For example, S = {|0, 1, +), |1, +, 0) , |+, 0, 1) , | — ,—,—)} in a three-qubit system, where 
we have defined |ib) = (|0) ib |l))/-v/2. It is straightforward to check that these basis states are 
orthogonal to one another. 

This set has the peculiar property that it is not possible to add an additional basis vector that 
is a product state. Suppose we add \a, b, c). This then requires that 

{a,b,c\0,l,+) = (a,6,c|l,+,0) = (a, 6, c|+, 0, 1) = {a,b,c\-,-,-) = 0. (F.l) 

But each of |a), and |c) can be, at most, orthogonal to any of the four states |0), |1), |+), and 
|— ). Hence, together, \a,b,c) can be orthogonal to three of the four states in S. This means that 
the set S cannot be extended, hence, the name unextendible product basis (UPB). Therefore, the 
subspace that is orthogonal to the space by the UPB contains no product states, nor mixture of 
them. 

What is the use of UPB? It turns out that it can be used to construct a bound entangled state. 
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Suppose S = {IV'i) , • • • J IV'n)} contains a UPB. Then the mixed state 

P=;o^ (|i-ElV'.)(V'.lj (F.2) 

is a bound entangled state, where D is the total dimension (e.g., D=8 in the above three-qubit 
example) . 

The state p is entangled because it lies in the subspace Ti. — 7is, in which, by construction, there 
is no product state, and hence, it cannot be written as a decomposition of pure product states. To 
see that its entanglement is bound, we can look at its partial transpose with respect to all bi-partite 
partitionings. As a product state is mapped to a product state under partial transpose and the 
identity is unchanged, p, under partial transpose, is mapped to a valid density matrix, which has 
non-negative eigenvalues. Therefore, p has PPT, and hence no entanglement can be established 
across any bi-partite cut via local operations and classical communication. We have then shown 
that the uniform mixture on the subspace complementary to that spanned by a UPB is a bound 
entangled state. In particular, the state 

/9 = ^(1-|0,1,+)(0,1,+| -|1,+,0)(1,+,0| -|+,0,1)(+,0,1| -|-,-,-)(-,-,-|) (F.3) 

is bound entangled. We have shown previously in Sec. 14.41 that such a state cannot violate, e.g., 
the Mermin-Klyshko-Bell inequality. 
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Appendix G 



Derivation of overlap of the ground 
state with the separable Ansatz state 



In this appendix we derive the results shown in Eqs. (j5.24p and (j5.25p . We first analyze h = 1/2 
case, viz., the even-fermion case. The lowest state |^i/2(^i^)) has zero number of Bogoliubov 
fermions. It is related to the state that has no c-fermions, i.e., = || • • • |) via 



|^i/2(r,/^)> = n cosC/'^(r,/i)e^*^-^-'^('-''^)^-''''^S^-^-i|J^) (G.la) 

m=0 

\^cos e jn{r, h) +i sin 9 rn{r,h)cl^c^j^^^^_^ \VL) . (G.lb) 



The Ansatz state is then 



m=0 



m)) = n[^o4 + ''''2 n (1-24q)(4-c.)] (G.2a) 

(G.2b) 

\n) (G.2c) 
^^■<'4'=? (G.2d) 

where we have suppressed the index (1/2). The term "^j^iC^jc] can be rewritten in momentum 



j=i i<i<j 



N 

= n( 


cos - + sm 










= cos^ 


2 




= cos^ 


— e 2 
2 


=ic]gtan2| 



space as 



— 2 — / 1\ 

E 44 = ^ E (G.3) 

l<j<l<N m=0 
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Thus, for even N 



j=l m=0 

whereas for odd N 



(e)) = (l + tan|^ct) n cos^|+-in^fcot^^^^cL4-_J|^^ (G.4a) 



|$(0)=(l + tan|^ct)cos| J] ^os^ | + ^ | cot ^^^^^44^1 1^^) • (G.4b) 

j=l m=0 V / 

Therefore, the overlap of the state \^i/2{r,h)) with fo^" even N is 



(*i/2(^,/i)l«>(e))= n I cos g^/^) (r, fe) cos^ I + sin ^Z^) (r, /^) sin^ | cot ""^"^^ , (G.5a) 



m=0 



whereas for odd 

(*i/2(r,/i)|$(0) =cos| n cos^^/2)(r,^) cos^ | + sin C^'H^, ^) sin^ | cot ^^^^^^ . 

m=0 \ / 

(G.5b) 

Next, we discuss the 6 = (odd-fermion) case. The lowest allowed state is the one with one 
'j^^ = c^^ fermion: 

|*o(r, h)) ^ |G(r, h)) = c^^^^ |G(r, h)) , (G.6) 
where \G{r, h)) is the state with no 7 fermions: 

|G(r,^)) = n [cos0io)(r,/i) +zsin0W(r,McWtcS^)tJ |Q) . (G.7) 

m=l 

Similar to the 6 = 1/2 case, by using 

l<j<l<N m=l 
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we obtain that for even TV 



m)) = (l + ViVtan i 4) cos2 | J] (^o^' | + \ ^c|„4_^ j , (G.9a) 

m=l ^ ' 

whereas for odd A'^ 

= (l + ViV tan i 4) cos | J] (cos' f + ^ sin^ | cot ^cLcjv^ ) |0) . (G.9b) 

m=l ^ 

Therefore, the overlap of |^'o(?", ^i)) with |4>(^)) for even N is 

(^'o(r-,/i)|$(0) = ViVsini n (cos4°)(r-,/i) cos^ | + sin ^ ^ /,) sin^icot^j , (G.IO) 
whereas for odd AT 



(^'o(r-,/i)|$(e)) = ViV 




(G.ll) 
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Appendix H 

Analysis of singular behavior of 
entanglement density 



In this Appendix we investigate the singular behavior of the entanglement density near the critical 
line h = 1 for both anisotropic (r 7^ 0) and isotropic (r = 0) cases. We begin with r 7^ case first. 

H.l Divergence of entanglement-derivative for the anisotropic 
XY models 



The starting point is Eq. (j5.26p . in which there is a maximization over the variable ^. The function 
to be maximized is 



F{i,r,h)= / d/i In [cos6l(/i,r,/i)cos2(^/2) +sin6'(^,r,/i)sin2(^/2)cot7r^] . (H.l) 



To find the stationarity condition, we demand the derivative with respect to ^ vanishes: 



d^F{i,r,h) 



2 Jq cos 6'(^, r, /i) 008^(^/2) + sin r, /i) sin^(^/2) cot vr/x 



cos r, h) — s'm9{fi, r, h) cot vr/i 



0. 



(H.2) 



Denote by (,*{h) the solution for fixed r. Then the field-derivative of the entanglement is 



dh£{r,h) = -±^dhF{e[h),h) = - A 



deih). 
dh 



d^F{^,h) +dhF{i\h) 



In 2 



dhF{e,h), 



(H.3) 



where the first term in the square bracket vanishes identically due to the condition ()H.2p . Thus 
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(dropping the * on ^ for convenience), 

dhSir,h) = -l^^dhFiCh) (H.4) 
In 2 



2 /■ 2 ^ dhCos9{fj.,r,h)cos'^{^/2) + dhsin6{fi,r,h)sm'^{S,/2)cotTTfi 
In 2 Jo cos6'(/i,r, /i)cos2(^/2) + sin0(//, r, /i)sin^(^/2)cot7r// 



Recall that tan 26{fj,, r,h) = r sin 2TTfi/{h — cos 27r^) and thus 



cos e = V(l + COS 26')/2, sin 9 = V(l^^cos2^j72, (H.6a) 

h — COS 2ttu , , 

cos20(^,r,/i) = . ^ .^"^ , (H.6b) 
y (r sm Ztt^)^ + (n — cos 27r/zj^ 



Putting everything in Eq. ()H.5p . we get 



^ ^ sin27r;U 
dh£{r,h) = -— / d/i' 



In 2 Jo (r sin 27r^)2 + (/i — cos 27r/i)2 



- {h - cos2tt fi) cos^{^/ 2) - \J ^/ + (/i - cos 27r//) sin^(^/2) cot vr/x 

/ ^ / = ) (H-7) 

y V + (/i — cos 27r^) cos2(^/2) + y V — {h — cos 27r/i) sin^(^/2) cot vr/x 



where V = (r sin 271^)"^ + {h — cos 27r//)2. 

We aim to explore the behavior near h = 1. First consider /i > 1 and define e = h — 1, which 
is the deviation from the critical point. Make the change of variables t = h — cos 2iti', giving lower 
and upper limits e and 2 + e, respectively. We further shift the integration variable by e, arriving 
at 

r /"^ 1 
dh£{r,h) = -^^J^ '^*(i_^2)i2 + 2(r2 + e)t + e2 

Vi^V - (t + e)cos2(e/2) - y^^^7"(t+^sin2(e/2)V2^ 

(H.8) 



ViyV +{t + e) cos2(e/2) + - (t + e) sin2(e/2)^/2^ 

where V = - r'^)t'^ + 2{r'^ + e)t + e^. 

As we inspect the limit /i — > 1 or e — > 0, we see that the above expression diverges, with the 
contribution coming from t small, i.e., infrared divergence. Large t{< 2) does not contribute to the 
divergence. Note further that only the second term in the numerator contributes to the divergence. 
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We then proceed to evaluate the integral by separating it into two parts: 



2 i-S p2 

= / + / , (H.9) 

Jo J 5 

with 5^1. In the first region we only need to keep t to first order at most. Also noting that for 
t, e <C 1, the first term in the denominator is much smaller than the second term, we get 



27r In 2 



-5 _i J^2(r2 + e)i + e2 + (t + e) 
dt 2 =■ (H.IO) 

2r^{t + ^) J^2{r^ + e)t + - {t + e) 



Next, we simplify the second term (ignoring e when there is no danger in doing so) 



V2(r2 + e)t + e^ + {t + e) Jt+ 



V2(r2 + e)t + e2 - {t + e) ^ 



^ + (H.11) 



Observing that only the first term on the right-hand side contributes to the divergence, we have 
that the divergent part is 



dt ^ = \ / dt^=^. (H.12) 



27rln2yo 2r2(t + ^) ^Tt irTTln2 Jq Vt+T Vi 

The divergence part is then (for 52r'^ /e'^ ^ 1) 



1 , / ^,^ \ 1 , ln(2\/^ 



2sinh-V'^2r2/e2 — In {2^/62r'^/€'^) = — lne + — ^ — (H.13) 



2^ 



4r7rln2 " ' 2r7rln2 V ^ 'J 2r7rln2 2r7rln2 

As the integral (|H.8P does not depend on the choice of 5, the part that involves 6 must be cancelled 
by the second half of the integration J^, which can be verified by direct evaluation. Therefore, for 
h very close to /ic = 1, we have 

d£ 1 

wr- —ln\h-l\. (H.14) 

9/i 27rrln2 ' ' ^ ^ 

In deriving the above divergence form, we have assumed that r ^ 0. Similar consideration can 
be applied to the case when h approaches 1 from below, and the behavior is the same. 
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H.2 Divergence of entanglement-derivative for the XX limit of 
the model 



We now analyze the r = isotropic case. It turns out that the analysis for this case is much simpler. 
To see this, we first note when r = we have the simplification 



cos 29{iJ,, h) = sgn(/i — cos 27r^). 



(H.15) 



The above expression changes sign when h = cos 27ri/. So let us introduce the variable = 
(cos^^ h)/{2iT). The expression for the entanglement density Eq. (|5.26|) becomes 

£{h) = —- — -max / d/i In [cos ^(/i, r, /i) cos^(^/2) + sin^(/x, r, /i) sin^(^/2) cot 7r//jH.16a) 
In 2 5 Jo 



max 

In 2 ? 



In sin (^/2) cot tt/U + / d/i Incos (^/2) 

J flQ 



Demanding stationarity with respect to ^ gives the condition 



(H.16b) 



[/"o - ^(1 -cos^)]sin^ = 0. 



(H.17) 



The solution ^ = gives the entanglement density for h > 1. It is straightforward to see from 
Eq. ()H.16aP that the entanglement density is identically zero. The solution fiQ — j{l — cos^) = 
gives 

cosC = 1- -cos^^/i. (H.18) 
vr 



This in turn gives the entanglement density for < h < 1: 



fio{h) In 



1 - 2fio{h) 



H — ln(l — 2/io(/i)) + / (i;U In cot vr^ 
2 Jo 



(H.19) 



Thus we have derived the entanglement density as a function of the magnetic field h in the XX 
limit. The result is shown in Fig. (j5.2p . 

We see from the figure that the entanglement density at /i = is the highest, which being 
£{0) = 1 - 27c/(7rln2) ^ 0.159 by evaluating Eq. (IHTO]) at /i = 0. The constant jc ~ 0.9160 
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is the Catalan constant. The entanglement density decreases monotonically as h increases until 
h = 1 beyond which it becomes zero identically. This qualitative behavior can be understood as 
follows. As the total z-component spin is conserved, increasing h simply increases the z-component 
spin of the ground state until h = 1 where all the spins are aligned with the field, hence there is no 
entanglement beyond this value of h. 

By directly taking the derivative with respect to h, we get 



dh£{h) 



1 



7rln2Vl - h'^ 



In 



cos ^ h 



n — COS" 



-1 



l + h 



hy 1-h 



(H.20) 



Near /i « 1, we have (putting 1 + /i = 2 and evaluating the limit in the argument of log function) 



dh ^ ' ' V27r 



{h 



(H.21) 



This completes our derivation of the singular behavior of the entanglement density near the critical 
points. 
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Appendix I 

Finite-size scaling 



The discussion here follows Barber [91]. In the vicinity of the bulk critical temperature Tc the 
behavior of a system should depend on y = L/^(T), where (,(T) is the bulk correlation length and 
L is the characteristic length of the system. How does the divergence of certain thermodynamic 
quantities emerge as the system size L grows? 

I.l Algebraic divergence 

Assume that some thermodynamic quantity at L — > oo diverges as t = (T — Tc)/Tc — > oo: 

Poo(T) ~ Coot"". (I.l) 

Finite-size scaling hypothesis asserts that for finite L and T near Tc, 

PUT) ~ rQp{l^'H), l^oo, i^O, (1.2) 

where / = L/a (a is some microscopic length), t = [T — Tc{L)]. The exponent oj can be determined 
by the requirement that the Pl(T) reproduces Poo{T) as / — > oo. Thus, 

Qp{x) ~ CooX~P, X ^ oo, (1.3) 

and io = p/v. We consider the case that the finite system does not exhibit a true transition, then 

Qp(x)^Qo, x^O. (1.4) 
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Prom this we have that at the peak or rounding temperature T^{1) (where Pl reaches the maximum 
or deviates significantly from the bulk value) 

Pi(r^(/)) ~ Qo/^/^ 1--00. (1.5) 

This means that the behavior of a thermodynamic quantity varies with the system size is determined 
by the bulk critical exponent. 

1.2 Logarithmic divergence 

Now assume the thermodynamic quantity P{T) diverges logarithmically as 

Poo(r) -Coolnt, t^O, (1.6) 

as in the field-derivative of the entanglement density for anisotropic XY spin chains. The finite-size 
scaling hypothesis in this case is to assume 

Pl{T) - Pl{To) ~ Qpil'/^i) - Qp(/V%), (1.7) 

where Tq is some non-critical temperature and to = {To — Tc{L))/Tc. The hypothesis has to recover 
the I — > oo limit at fixed T, which requires 

Qp{x) ^ Coo In 3;, X —>■ oo. (1.8) 

Thus in the limit f — >^ at fixed large Z, we have 

Pi(rc(L))~-^ln/ + 0(l), (1.9) 

\iQp{x) = 0(1) as a; — ^ 0. This allows us to obtain the exponent v by analyzing how the divergence 
develops as the system size I (which is in our spin-chain entanglement problem) increases. 
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